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Abstract. We consider the stationary solutions for a class of Schrodingcr 
equations with a symmetric double-well potential and a nonlinear perturba- 
tion. Here, in the semiclassical limit we prove that the reduction to a finite- 
mode approximation give the stationary solutions, up to an exponentially small 
term, and that symmetry-breaking bifurcation occurs at a given value for the 
strength of the nonlinear term. The kind of bifurcation picture only depends 
on the non-linearity power. We then discuss the stability /instability proper- 
ties of each branch of the stationary solutions. Finally, we consider an explicit 
one-dimensional toy model where the double well potential is given by means 
of a couple of attractive Dirac's delta pointwise interactions. 



1. Introduction 

Here, we consider the stationary solutions of the nonlinear Schrodinger (hereafter 
NLS) equations 

ihj£=H Q 4, + eg{x)\1>\*'1>, = 1, (1) 

2 



where e e R and || • || denotes the L 2 norm, 

flb— AjA + V, A = S * (2, 

is the linear Hamiltonian and g(x)\ip\ is a nonlinear perturbation. For the sake 
of definiteness we assume the units such that 2m = 1. 

Atomic Bose-Einstein condensates (BECs) are described by means of nonlinear 
Schrodinger equations of the type (JTJ) where Hq represents the Hamiltonian of a 
single trapped atom and the nonlinear term |f/>| 2cr , a = 1, 2, . . ., is the (a + l)-body 
contact potential [35]. In fact, BECs strongly depend by interatomic forces and 
the binary coupling term |'0| 2 ' ! /' usually represents the dominant nonlinear term and 
equation |T]) takes the form of the well-known Gross-Pitaevskii equation [35] . Even 
if in most of the applications the parameter a takes only integer and positive values, 
here we take that a can assume non integer values too, as considered in [33 . It is 
worth mentioning also the fact that equation (JlJ with nonlinearity corresponding 
to the power-law |V , | 2<T ! where the parameter a takes any positive real value, is used 
in other contexts, including semiconductors [26 and nonlinear optics 6, 34, 35 . 



Date: January 19, 2013. 

One of us (A.S.) is very grateful to Giuseppe Mazzuoccolo for useful discussions about the 
Budan-Fourier theorem, and to Riccardo Adami and Hynek Kovarik for useful discussions on NLS 
equations with singular pointwise interactions. 

1 



2 



REIKA FUKUIZUMI AND ANDREA SACCHETTI 



In this paper we consider the case of symmetric potentials V with double well 
shape; the function g(x) is a bounded regular function (in the following we assume, 
for argument's sake, that g(x) has the same symmetric properties as V(x)). 

If the nonlinear term is absent then the linear Hamiltonian Hq has even-parity 
and odd-parity eigenstates: the d-dimensional linear Schrodinger equation with a 
symmetric double well potential has stationary states of a definite even ip + and 
odd-parity <ys_ , with associate nondegenerate eigenvalues A+ < A_ . 

However, the introduction of a nonlinear term, which usually models in quantum 
mechanics an interacting many-particle system, may give rise to asymmetrical states 
related to spontaneous symmetry breaking phenomenon. 

In the semiclassical limit and in the two-level approximation has been seen |32j 
that the symmetric/ antisymmetric stable stationary state bifurcates when the adi- 
mensional nonlinear parameter r/ takes absolute value equal to the critical value 

rf = 2"/a. (3) 

The parameter rj is associated with the coupling factor of the nonlinear perturbation 

by 

rj = ce/oj (4) 

and it is the effective nonlinear coupling factor, where u is the (half of the) splitting 
between the two levels 

W =~(A_-A+) (5) 

and c is a constant defined below in §2.2. In fact, in the semiclassical limit (or also 
for large distance between the two wells) the splitting uj is exponentially small, as 
h goes to zero. Furthermore, in [53] it has been also seen that for a less than a 
critical value 

_ 1 

^threshold ~Z 

then a supercritical pitchfork bifurcation occurs; on the other hand, for a bigger 
than the critical value cr t hreshoid a subcritical pitchfork bifurcation associated to 
the appearance on a couple of saddle node points occurs. 

It is worth mentioning the fact that the main problem consists in proving the 
stability of the two-level approximation (which basically is a two-mode problem) 
with respect to the NLS equation ([!]). So far, the stability of the two-level approxi- 
mation has been proved, in the semiclassical limit, only for times of the order of the 
beating period T = 2irh/uj [51] . or for exponentially large times (that is of the order 
e T ) under further assumptions as proved by [5J. In fact, our previous approach 
was rather efficient in order to study the dynamics, but only give a partial result in 
order to look for the stationary solutions. Recently, Kirr, Kevrekidis, Shlizerman 
and Weinstein [53] has considered the stationary solution problem for the Cauchy 
problem (fT]) with h fixed (i.e. H = 1) in the limit of large barrier between the two 
wells, and in the case of cubic nonlinearities. In their seminal paper they make use 
of the Lyapunov- Schmidt reduction method to the two-level approximation equa- 
tion for the stationary solutions. In such a way they overcome the limit of the 
method applied by [51] for the study of the stationary solutions. Furthermore, 
they also applied the same method in order to study the orbital stability of the 
obtained solutions. 
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In this paper we follow the ideas developed by [23] , adapted to the semiclassical 
limit and considering the case of any positive and real nonlinearity power ct, in order 
to study the stationary solutions of equation ([T]) and their stability properties as 
function of the nonlinearity power a. In particular we are able to prove that the 
result obtained by |32j for the two-level approximation, concerning the existence on 
the critical value (^threshold, holds true for the whole Cauchy problem ([1]), too. To 
this end we prove the stability of the two-level approximation, when restricted to 
the stationary problem, and then we apply a generalization of the Budan-Fourier 
theorem [TU] in order to count all the branches associated to the stationary solutions. 

It is worth to mention the fact that the stability of the two-level approximation 
holds true in order to classify the stability/instability properties of the stationary 
solutions, too. In fact, stability/instability properties of the stationary solutions 
for the two level approximation are easily obtained since such an approximation 
has a finite-dimensinal Hamiltonian structure. On the other side, orbital stabil- 
ity/instability properties of the stationary solutions of the full nonlinear problem 
are much harder to obtain. However, in this paper, by making use of the meth- 
ods developed by Grillakis, Shatah and Strauss [El US], and succesfully applied by 
[23] for double well problems with cubic nonlinearity, we prove the equivalence be- 
tween the stability/instability properties when we restrict our problem to the case 
of attractive nonlinearity and when we restrict our analysis to the " ground state" . 

There are already many studies on the existence of stationary solutions and 
the stability of Eq.([T]) in the semiclassical limit (e.g., [13l [T8J HQ])- However, our 
aim is to understand what happens with double-well problem. When we con- 
sider the stationary problem with symmetric double- well and nonlinearity strength 
large enough, the bifurcation picture tells us that we have asymmetrical stationary 
solutions localized on just one well, as well as asymmetrical stationary solution 
delocalized between the two wells. The first type of solution was obtained, but 
the second type of solution was not considered in [13 , and it is identified with the 
multi-bump stationary solution studied in, e.g., Also it would be important to 

understand the destruction of the beating motion in the framework of the dynamics 
(see [17] for related topics). 

The paper is organized as follows. In Section[2]we recall some preliminary spec- 
tral results for Schrodinger operator with double well potential in the semiclassical 
limit, we introduce the main assumptions and we collect some general global well- 
posedness results for the Cauchy problem ([1]) . In Section [3] we prove (Theorem 
[TJ concerning the occurrence and the nature of spontaneous symmetry breaking 
phenomenon for equation ([1]) by applying, in the semiclassical limit, the Lyapunov- 
Schmidt reduction method to the two-level approximation and some results of the 
theory of numbers in order to count the number of solutions of a polynomial-type 
equation coming from the two-level approximation. In Section |4] we consider the 
dynamical properties of the stationary solutions of the two-level approximation, 
which has Hamiltonian form. In Section [5] we consider the orbital stability proper- 
ties of the ground state stationary solutions. Appendix is devoted to an application 
of all the arguments in the previous sections to an explicit one dimensional toy 
model where the double well potential is given by a couple of attractive Dirac's 
delta interactions. 

Notations. Hereafter, 
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• y = 0(x), means that for any < a < 1 there exists a positive constant 
C := C a such that \y\ < C a \x\. Here, as usual y = 0(x) means that there 
exists a positive constant C such that \y\ < C\x\, and x ~ y means that 
lim^o | = C for some Cel; 

• || • ||p and || • || denote the norm of the spaces LP and L 2 , (4>,<p) = J 4><p 
denotes the scalar product in the Hilbert space L 2 ; 

• C denotes any positive constant which value is independent of %. 

2. Main assumptions and preliminary results 

Here, we recall some preliminary results. Throughout the paper we always 
assume the Hypotheses below in this section. 

2.1. Linear operator. Here, we introduce the assumptions on the double-well 
potential V and we collect some well known results on the linear operator Hq. 

Hypothesis 1. The potential V(x) is a bounded real valued function such that: 

i. V is a symmetric potential. For the sake of definiteness we can always as- 
sume that, by means of a suitable choice of the coordinates, V is symmetric 
with respect to the spatial coordinate x\, that is 

[S,V]=0 (6) 

where 

[Sip] (xi,x 2 , ...,x d ) = 1p(-Xi,X 2 , ■ ■ .,x d )- 
Hence, the Hamiltonian Hq is invariant under the space inversion: [S, Ho] = 
0, 

ii. V € C°°(R d ); 

iii. V(x) admits two minima at x = x±, where X- = Sx+ =/= x + , such that 

V(x)>V mm = V(x ± ), Vx£R d ,x^x ± . (7) 
For the sake of simplicity, we assume also that 

VV{x ± ) = and Hess V(x±) > 0. 

iv. Finally we assume that the two minima are not degenerate: 

V x = Umirjf| a .|_ >0O V'(a;) > V min . (8) 

Remark 1. In fact, some assumptions on V may be weakened. In particular, 
the case of degenerate minima, that is det [Hess V(x±)] — 0, could be treated in a 
similar way; however, we don't dwell here on such details. Furthermore, bound- 
edness of V is assumed just for sake of definiteness: if V is not bounded we could 
make use of the argument by [3] in order to prove the well-posedness of the Cauchy 
problem under some assumptions of the behavior of the potential at infinity. 
For instance, we could assume that there exists a positive constant < m < 2 such 
that for large \x\ 

C(x) m < V{x) < C- X {x) m , (x) = (1 + \x\ 2 ) 1/2 , 

for some C > 0, and 

d 

\d% . ..d^V(x)\ < C a (x) m -M , \a\ =Y,<*i> 

3=1 
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for any multi-index a £ N d . 

The operator H formally defined by © admits a self-adjoint realization (still 
denoted by Ho) on H 2 (R d ) since V is a bounded potential. 

Let <j(Hq) = adUa ess be the spectrum of the self-adjoint operator Hq, where ad 
denotes the discrete spectrum and o~ ess denotes the essential spectrum. It follows 
that 

o~d C {V min ,V~) and a ess = [V~, +00) . 

Furthermore, for any h £ (0, H*), for some h* > fixed and small enough, it 
follows that ad is not empty and, in particular, it contains two eigenvalues at least 
Ai_ and Al where \\_ < A 1 , and 

inf [C-Xi]>CH, (9) 
for some positive constant C independent of ft. 

Remark 2. Actually, from Hypothesis^ and for h small enough in general it follows 
that for some E > V m in then 

Od n (V min ,E) 

is given by a sequence of couple of nondegenerate eigenvalues \±, j — l,2,...,n 
where n ~ hr 1 , such that X 3 , < AiL and 



inf 



c-Ai 



> en (10) 



hold true. In fact, degeneracy may occur for some j > 1 only in special cases, for 
instance when other symmetry properties for the potential V are present (see, e.g., 
[20j ). Hereafter, for the sake of definiteness, we assume that degeneracy does not 
occur and that \1U\) holds true for any j — 1, 2, . . . , n. 

Let tp± be the normalized eigenvectors associated to A^_ , then (p± can be chosen 
to be real-valued functions such that 

Sifi = ±<4; (11) 

Furthermore 

Lemma 1. The eigenvectors (p± belong to the space H 2 (M. d ) n L p (M. d ) where 

{< +00 ifd=l 
< +00 ifd = 2 . (12) 

<2d/(d-2) ifd>2 

In particular, it follows that 

||Vv4ll < Cjtr 1 ' 2 and ||v4llff 2 < C i h ~ l (I 3 ) 

and 

W^Wp < C 3 h- d2 ^ , (14) 
for some positive constant Cj, independent on H. 
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Proof. Indeed, Lp± is normalized and it satisfies to the following eigenvalue equation 
— h 2 Aip± — (A± — V")y±, from which immediately follows that 

((Ai-FV ± ,4> L2( oi) 



< 



2 



where 



n 3 ± ={xeM d \ V{x) < A J ± } 



is such that A± — V > \± — V m in > Cjh for any hxed j and H small enough. 
Similarly 

ftW±ii 2 = |(Ai- v)cpi\\ <<w±n. 

since V is a bounded potential. Estimate (fl4|) follows by means of the Gagliardo- 
Nirenberg inequality: 

where S = ^H^d. □ 

Remark 3. Actually, tf± € LP for any p and, by means of the Riesz-Thorin 
interpolation Theorem, inequality holds true for any p independently on the 

dimension d (see, e.g., [31 J. Indeed, by means of the semiclassical expression of 
ipj it follows that || v^j ||oo < Cjh" d ^ 4 . 

The splitting between the two eigenvalues 

J = |(Ai - A J + ) (15) 

vanishes as h goes to zero. In order to give a precise estimate of the splitting uji 
we make use of the fact that V is a symmetric double- well potential with non-zero 
barrier between the wells. That is, let j be fixed and let 



inf 

7 



y/[V(x)-V min ] + dx > 0, (16) 



be the Agmon distance between the two wells; where 7 is any path connecting the 
two wells, that is 7 € AC([0, 1], R d ) such that 7(0) = x_ and 7(1) = x+, and where 
[•]_!_ = max(-,0). From standard WKB arguments (see [20J for details) then it 
follows that the splitting is exponentially small, that is 

^ = b{e- p / h ) . (17) 

Let LfP R L be the normalized single well states associated to the linear eigenstates 
(p± by means of 

<? R = {<P i + + <pL)/y/2 (18) 

and 

^ = (^ + -^l)/A (19) 
They are localized on one well in the sense that and for any p £E [2, +00] then 

\\^ L \\p=0(e-"/ h ). (20) 



BIFURCATION AND STABILITY FOR DOUBLE-WELL NLS 



7 



More precisely, these functions are localized on only one of the two wells in the 
sense that for any r > there exists c := c(r) > such that 

f \^ R {x)\ 2 dx=l+0{e- c ' h ) 

JD r (x + ) 

and 

\p> L (x)\ 2 dx = 1+ 0{e- c / h ) 

D r (x— ) 

where D r (x±) is the ball with center x± and radius r. For such a reason we call 
them single-well (normalized) states. 

Remark 4. In the following, for the sake of definiteness, we restrict ourselves to the 
couple of eigenvalues \\_ and X_, corresponding to the lowest energies. Hereafter, 
we simply denote them by A± dropping out the index 1, and tp± denote the associated 
eigenvectors. The symmetric solution (p+ is the first eigenfunction of Hq, so it is 
positive. We remark that the existence of the stationary solutions for {!]) and their 
dynamical stability still hold true when we consider all the unperturbed energy levels 
A± provided that degeneracy does not occur as discussed in Remark^ 

2.2. Assumption on the non linear term. In order to obtain some a priori 
estimates of the wavefunction \ip\ 2cr i^ we introduce the following assumption on the 
nonlinearity power a. 

Hypothesis 2. We assume that 

n ^ „ ^ / +°° if d=l, 2, , . 

i — - — if d > 2 ( 21 ) 



i-2 



where d is the spatial dimension. 
Let 

C R = {^+\g^ x ) and C L = {^+\g^+ l ) 

where Cr — Cl because of the symmetric properties of g and V. We assume also 
the following scaling limit. 

Hypothesis 3. Letui = i(A_ — A + ) be the splitting U5)) satisfying to the asymptotic 
estimate |J7| ). We assume that the real-valued parameter e depends on h in such 
a way 

\rj\<C where 7} = — , c:=C R = C L , (22) 

for some positive constant C , independent of h. The parameter rj plays the role 
of effective nonlinearity parameter. Hereafter, we assume that g{x) has the same 
symmetry property @) of the potential V and it is such that (ip R +1 , g(p R +1 ) ^ 0. 
In particular, for the sake of definiteness, let 

K + \<?^ +1 )>0. (23) 

2.3. Existence results in H 1 and conservation laws. The results below follow 
from [5] and from the a priori estimate given by [31] , 
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2.3.1. Local existence in H 1 . Let the initial state ip Q E H 1 , then there exists T* > 
and an unique solution ip(x,t) € C([0, T*), H 1 ) n C 1 ([0, T*), H' 1 ) of Q}, where 
T* = +oo or IIV^H — > +oo as t — >• T* — 0. Furthermore, the conservation of the 
norm and of the energy hold true for t £ [0, T*]: 

U(;t)\\ = U°(-)\\ 

and 
where 

cr + 1 

represents the energy functional. 

2.3.2. Global existence. The solution ?/> of (JTJ globally exists, that is T* = +oo, 
provided that the state is initially prepared on the first N states of the linear 
problem, for any TV fixed, and h is small enough. Indeed, this fact immediately 
follows from a priori estimate of the norm of the gradient of the wavefunction [31] . 

Remark 5. The solution ij}(x, t) globally exists for both positive and negative values 
of the parameter e, provided that H is small enough and e satisfies Hyp. That 
is, because of the scaling assumptions, blow-up effect cannot occur. 

3. Stationary solutions and bifurcation 

Since the beating period T — 2^5 plays the role of the unit of time it is convenient 
to introduce the slow time 

Ljt 

T = Y' 

then equation ((TJ takes the form (here ' denotes the derivative with respect to r 
and where, with abuse of notation, -0 = tJj(t,x)) 

i^ 1 ^ H ^ + eg\^\ 2 ^ , ||V>G,t)|| = 1. (24) 

In order to study the stationary solution we set 

i/>(x, t) = e~ lXT ^4,{x), = 1 , A = Q + uE, 

where 

n = \ [a+ + a_] . 

As specified in Remark S] we restrict ourselves, for the sake of definiteness, to the 
first couple of energy level A^j_ , where we simply denote them by A± dropping out the 
index 1; similarly ip± denote the associated eigenvectors and v?-R,£ the associated 
single- well states. 

Hence, equation (|24l) takes the form 

Xij = H Q iP + eg\ib\ 2 °iP, ||V(-)II=1- (25) 

Now, let us set 

ip(x) = a R ip R (x) + a L cp L (x) + ip c (x) , (26) 

where 



ip c (x) = U c ip(x) 
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and 

dR = ((pR,tp) and a L = (tp L ,il;) 
are unknown complex- valued values. Here, 

n c = i-n, n = [^ + , .) v+ + (<p_,.)<p_] 

denotes the projection operator onto the eigenspace orthogonal to the bi-dimensional 
space associated to the doublet {A±}. 
Since 

H ip = aRH ifR + ciLHoipL + H ip c 

= a R [VtipR - ujip L ] + a L [-ucpR + tttp L } + H Q ip c (27) 

then, by substituting (l26l) in (|25[) and projecting the resulting equation onto the 
one-dimensional spaces spanned by the single-well states ifR and ifL, and on the 
space U c L 2 (R d ) it follows that (|2"5]| takes the form 

Ea R = -a L + r R , r R = r R {a R ,a L ,^ c ) = ^(ip R , g\tp\ 2,7 ip) 

Ea L = -a R + r L , r L = r L (a R , a L ,ijj c ) = ~{<fL, g\^\ 2a ^) (28) 

£Vc = h i H o ~ fi]V»c + r c , r c = r c (a R ,a L ,ijj c ) = ^n c5 |V| 2<7 V 
with the normalization condition 

\a R \ 2 + \a L \ 2 + \\M 2 = 1- 

Remark 6. Since equations 128\) has stationary solutions \26\) defined up to a phase 
term then we can always assume, for the sake of definiteness, that the stationary 
solution of equation H25\) is such that is a real-valued positive constant: or G C 
and ol G K + . Furthermore, we remark that [Hq,S] = and [g,S] — 0; hence, if 
ip is a stationary solution of equation (25\) associated to a given value X, then Sip 
is a solution associated to the same level, too. 

Then, collecting the results from Lemmata 3 and 4 (and the associated remarks) 
by [31] we have the following. 

Lemma 2. Let p be the Agmon distance between the two wells defined as in I116\) . 
It follows that 

rR,L(aR,a L ,ip c ) = r RiL (a R ,a L ,0) + r RL (a R , a L , ip c ) 

where 
(i) 

rR,L(aR,a L ,0) = -C R , L \a R , L \ 2 °a R , L + ~0(e-P' h ) (29) 

and 

Cr, l = { V R,L,gW R M 2a VRL) = ^tl9^l) = o(h- d ^ 2 ) ] (30) 

by the symmetry assumptions it turn out that 

Cr — Cl ■ 
(ii) The remainder terms are estimated as follow 

\i% L \ < -ch-www 
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where 



7 i 1 + (2 - d)7 i/ d > 2 • ^ J 



Here we come with the existence result of stationary states for the nonlinear 
Schrodinger equation (|2"5j) . 

Theorem 1. Let 

afl = pe l9 , aL = q and z = p 2 — q 2 (32) 

where p,q € [0, 1] and € [0, 2ii). Let h G (0, h*), where h* is small enough, let p 
be the Agmon distance between the two wells and let r\ be the effective nonlinearity 
defined by H22]) . Then the stationary problem (28\) always has 

- a symmetric solution ip s E such that 

9 s = d{e-P^/ h ), z s = 0{e~ pl/h ), 

associated to 

- an antisymmetric solution ip^ such that 

9 a =ir + 0{e-^' h ), z a = d(e- p ^ h ), 
associated to 

E:=+l + V ^ + d(e-P^ h ). 

Furthermore, in the case of negative (resp. positive) rj, then asymmetrical so- 
lution ipE* corresponding to 9 as — 0{e~ p ~</ h ) (resp. 9 as = n + 0{e~ pl / h )) may 
appear as a result of spontaneous symmetry bifurcation phenomenon. That is: 

- for a < o-fhreshoid the symmetric (resp. antisymmetric) state corresponding 
to z s — 0(e~ P7 /' i ) bifurcates showing a pitchfork bifurcation when the adi- 
mensional nonlinear parameter \rj\ is larger than the critical value rf given 
by (see Fig. [7J panel (a)) 

V = — 

a 

- for a > crthreshoid two couples of new asymmetrical stationary states 
appear as saddle-node bifurcations when \r]\ is equal to a given value rf~ 
such that 77 + < rf ; then, for increasing values of \rf\ two branches of the 
solutions disappear at \r}\ = rf showing a subcritical pitchfork bifurcation 
(see Fig. [7J panel (b)). The critical value r/ + is given by rj(z + ) where 

(33) 



l + z\° \ 



( \ 2z 

v[z) = VT=z- 

and z + € (0, 1) is the non zero solution of the equation rj'(z) = 0. 
In all the cases, the remainder term ip c of the stationary solutions is such that 

WcWh* = 0{e- p ' h ) . (34) 
The critical value crthreshoid is given by 

_ l r 

^threshold ^ 
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and it is an universal value in the sense that it does not depend on the shape of the 
double well potential as well as on the dimension d. 

Remark 7. Concerning the symmetric solution ip E = a^ipn + a^ipL+ipc we should 
remark the above statement implies that the corresponding level E is non- degenerate 
in the sense that we have only this stationary solution corresponding to such value 
of E. On the other side, by means of a symmetrical argument, then Sip s E = 
a RfL + 0'L i PR +$?Pc is a solution associated to same level E, too. Hence, ip E and 
Sip s E coincide, up to a phase factor. From this fact and from Theorem [JJ it turns 
out that 9 s and z s are exactly zero: 

9 s = and z s = . 

Similarly, it follows that 

9 a = 7T and z a = 

and 

9 as = (respectively 9 as = tt) 

for negative value of r\ (resp. for positive value of rj). By means of a similar 
argument we can also conclude that the stationary solution is, up to a phase term, 
a real valued function; indeed if tp is a solution associated to a given level E, then 
ip is a solution associated to the same value E, too. 

Remark 8. Because of the technical assumptions on a, this critical value crthreshoid 
makes sense for the non-linear Schrddinger equation 125\) only in dimensions 1 and 
2. This is not the case when we restrict our analysis to the two-level approximation. 

Remark 9. From Theorem\J\it appears that we have only two pictures, accordingly 
with the value of a. In Fig. [JJ (panel (a)) we consider the bifurcation scenario for 
the imbalance function z — \aj^\ 2 — \a^\ 2 appearing when a < a 'threshold- In Fig. 
[7] (panel (b)) we consider the bifurcation scenario appearing when a > o~ threshold • 
The same picture has been previously obtained for the two-level approximation (see, 
e.g., [32) ) where we have taken ip c = 0; in fact, ip c is exponentially small as proved 
in Theorem]]} 

Remark 10. The stationary solutions tp := ipE, associated to the level E, given in 
Theorem [JJ are such that 

||V^M| < CVA (35) 

and 

\\iPe\\p < CA dE ^ (36) 
where p satisfies condition h!2\) and where 

H(lp E ) ~ V m in fe _l 



A = 



ft 2 



is the energy functional defined on H 1 (K d ) n L 2( - a+1 1 (R d ). Indeed, estimates fffSj) 
and \30j) hold true for any vector tp belonging to the space II(L 2 ) (see Theorem 2 
in |31j ). The results finally follow from this fact and since H c ipE = 0(e~ p / h ). 
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Figure 1 . In this figure we plot the graph of the stationary states 
of the non-linear Schrodinger equation (|25[) as function of the non- 
linearity parameter 77 for nonlinearity a = 1 < athreshoid (panel 
(a)) and for nonlinearity a = 5 > <J threshold (panel (b)); here 
z = \o-r\ 2 — l a i| 2 is the imbalance function. Full lines repre- 
sent stable stationary states and broken lines represent unstable 
stationary states, where the notion of stability is referred to the 
dynamical stability associated to the Hamiltonian system given by 
the two-level approximation, as discussed in S}3J and also to orbital 
stability, as discussed in ^5] in the case of attractive nonlinear case 
(i.e. 77 < 0). 

3.1. Proof of Theorem[TJ Here, we prove the existence of the stationary solutions 
by making use of the Lyapunov-Schmidt method and applying some results of the 
theory of numbers in order to count the number of stationary solutions of the 
equation coming from the two-level approximation. In this section, for argument's 
sake, we take 77 > 0; however, the same results still hold true also for ?y < 0. 

Lemma 3. We consider the following equation 

[H -n-wE]ip c + eIL c g\ip\ 2a il} = Q. (37) 

where the nonlinearity power a satisfies condition \21\) . For any fixed C > let 

D = {{a Rl a Ll E) e C 2 x M : \a R \ 2 + \a L \ 2 < 1, \uiE\ < Ch 2 } . 

Then, for there exists H* > small enough such that for any fig (0, K*) then there 

exists an unique solution ip c G H 2 of equation J^T] ) depending on a R , and E, 
and such that 

max HcWm = (e- p/h ) , as % -> . (38) 

(a R ,a L ,E)£D \ / 

Proof. Recalling that 

ip = <p + ip c , where, ip = a R ip R + a L ip L , 
then (|3"?| takes the form 



= F{ip c ) 



(39) 
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where 

F{$ c ) := F(ip c ; a R , a L ,E) = -e [flo - « - uE]' 1 n c .g|Vf CT ^ (40) 
and where 



\H -n-uE] n c <Cift" 1 (41) 

£(£2-KfP) 

for some positive constant Ci and for ft small enough, since and since = 
0(h 2 ). On the other side we have that 



\\F{u)-F{v)\\ H2 < e ^\\\f\ 2 y-\g\^g\\ 



,^2 in j.j2<r |„|2er, 

< ^iKi/r+| 9 m/- 9 |ii 

< ^(11/11^+11^011/-^ 

for some positive constant C2, where we set 

f = ip + u and 5 = ip + u, 

with \clr\ 2 + |ol| 2 + ||u|| 2 = 1, |ai?| 2 + |aL| 2 + ||w|| 2 = 1. We have indeed made use of 
the Holder inequality and of the Gagliardo-Nirenberg inequality with a satisfying 
condition ([2~Tj): if 2pa < b and 2p/(p — 2) < b where b — +00 if d = 1,2 and 
b = 2dj (d — 2) if <i > 2, i.e. a satisfies (|2~T1) . Finally, we get the wanted estimate 

\\F(u) - F(v)\\ H2 < e ^^{ m ax[||^ + u || H2 ,||^ + «|| ff2 ]} 2CT Hu-wll^ (42) 

provided that cr satisfies condition (f^Tj) . 
Now, let C 3 = max[Ci,2 2,T C 2 ] and let 



f 


ft 


l/2<r > 




2 2 -+ 2 3C 3 e_ 


>IM|i?a | 



K = {u e H 2 : \\u\\ H 2<c{h)), c(ft) = max- 

r "i 1 / 2<T 

Since \\<p\\h 2 = 0(h~ 1 ), by LemmaQ] and e = 0{e~ p / h ) then c(ft) = 2 ^+^ 3C - 
Then .F is an operator from K to if; indeed, from (|4T|) and (|42|) it follows that 

ll^WlUa < e C 3 ft- 1 ||u + ^|| 2 7+ 1 < [2eC 3 ft- 1 (2c) 2CT ] c(ft) = ic(ft) < c(ft) . 

Moreover, -F(u) is a contraction in K: 

\\F(u) - F(v)\\ H 2 < CzetT 1 [2c(K)] 2a \\u - v\\ H * <\\\u- v\\ H * ■ 

Hence, equation 

F(u) = u 

admits a unique solution ip c in K for any (aij, a^, E) £ D and any e satisfying Hyp. 
[3] This solution is given by the limit of the following sequence {u„}^L where 

u = and u n+1 = F(u n ) . 

In particular (the convergence is in H 2 ) 

+00 +00 
i'c= lim u n = ^ [uj+x - uj] = V [F(uj) - F(uj-i)] 

3=1 3=1 
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Since 



\\F(u j+1 ) - F( Uj )\\ H2 < Ckeh- 1 [2c(h)] 2 ' T \\F(u j )-F(u j - 1 )\\ H2 



< 



[C 3 eft- 1 [2c(/i) 



\F(uo)\\ 



H 2 



then wc have that 

Wc\\h* < 



\F(*o)\\ 



H 2 



< 



I 2(7+1 



1 - C 2 eh- 1 [2c(h)} 2 ° 

= 0{e- p ' h ) (43) 
Since the constants C\ and Ci depend on a R , ol and E in such a way that 



and 



max C\ < +oo 

\u)E\<Ch 2 



max C*2 < +oo 

\a R \ 2 + \a L \ 2 <l 



then the estimate (|43l) uniformly holds true on the set -D. 



□ 

Remark 11. By means of the same arguments it follows that ip c G H , as function 
on an, ol and E, admits the first derivatives and in particular these derivatives 
satisfy estimate (38\) in the sense that 



max 

(a R ,a L ,E)eD 





dipc 




dip c 




dip c 






dE 


H 2 


da R 


7 

H 2 


da L 


H 2 - 



= o 



(e- p/r ^j , as H^O. (44) 



We can also give an estimate of the dependence of ip c on the parameter e; this 
estimate will be given in Lemma^ 

Now, setting ip c = ij) c 'a R , cll,E) in ([28]). let any < p' < p fixed, let 

v = e r p 'il n (45) 

where 7 is defined in equation ([31]) , and making use of Lemma [2] then ([28]) takes 
the form 

Ea R = -a L + ri\a R \ 2cr a R +vf R {a R , a Ll E) 
Ea L = -a R + r]\a L \ 2,7 a L +vf L (a R ,a L ,E) (46) 
1 = |afl| 2 + M 2 +vf c (a Rl a L ,E) 

where f R , f^ and f c are uniformly bounded on D with their first derivatives. Since 
Lemmalland Remark [TTl and recalling that e/uj = r/ / (ip R +1 , gip R +1 ) = 0(hr da l 2 ). 
From then gB]) takes the form 

Ep = -qe- i9 + rip 2a+1 +ve- ll> f R 
Eq = -pe ie + riq 2a+1 +uf L 
1 = p 2 + q 2 +vf c 

By taking the real and imaginary part of the previous equations we obtain the 
following system 

G(p,q,E,8;v)=0 (47) 



D' = {(p,q,E,6) e [0,1] 2 x E x [0, 2tt) : p 2 + q 2 <l, \wE\ < Ch 2 } 
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and where G — {G\, G2, G3, G4) are given by 
1 



Gi 



E 



[-2pqcos9 + r,(p 2a+2 + q 2 °+ 2 ) + ^(pe^ 6 ' f R + qf L )] 



= E + 2pqcos9 - i](p 2a+2 + q~ 
G 2 = (p 2 + q 2 ) sin 9 + u^{e- ie f R - pf L ) = (p 2 + q 2 ) sin 9 + vf 2 
G3 - (p 2 -q 2 )cos9 + r ] pq(p 2 ° -q 2 °) + v$t(qe- ie f R -pf L ) 
= (p 2 -q 2 )cos9 + T)pq(ip 2 ° - q 2 °) + i// 3 
q 2 -I- I' J, - 1 = p -r '/ - 1 -i- i') 



G 4 = p 2 +q 2 + vf c -l=p 2 +q 2 - l + vfi 



where fj, j — 1, 2, 3, 4, are uniformly bounded on the set D' with their hrst deriva- 
tives. 

From equations Gi = and G4 = we obtain that 

p 2 + q 2 = 1 + 0(v) and 9 = 0{v) , 9 = n + 0{v) . 

From this fact and from equations Gi = and G3 = we finally obtain the 
equations 



G± + 0(v) =0 



2<j+2 



,2cr+2 



(48) 
(49) 



where the asymptotics is uniformly on D' , the index + corresponds to the choice 
9 = 0{v), the index — corresponds to the choice 9 = ir + 0(v) and where 

G ± = ±[(p 2 -q 2 )±r)pq(p 2 °-q 2 °)] ■ 
The imbalance function z = p 2 — q 2 is such that 



P 



l + z 



0(v) and q 



1 - z 



2 y ' * V 2 

and thus equations (|48|) and (|49| take the form 

f±(z,v) + O(v) = 

E± = TV^~z- 

where 

" 'l + z 



O(iz) 



<T + 1 



O(iz) 



f±{z,v) =z±r) 



(50) 
(51) 

(52) 



Since the asymptotic term 0{v) in (|50p . with its derivative with respect to z, is 
uniform with respect to z G [—1, +1] then it is enough to look for the solutions of 
equations f±(z,rj) = 0. 

Of course, equation 

f±(o,v) = o 

holds true for any 77; that is the symmetric stationary solution (z = 0, 9 = 0) which 
is positive and the antisymmetric stationary solution (z = 0,9 = tt) exist for the 
nonlinear problem (up to an exponentially small perturbation) as well as for the 
linear one. 



1(5 



REIKA FUKUIZUMI AND ANDREA SACCHETTI 



Since we have assumed, for the sake of definiteness, 77 > 0; then equation 
f+(z,T)) — does not have non zero solutions, indeed the derivative of /+ with 
respect to z is given by 



1 + z \ 1 1 — z 



f'4z,v) = 2 [1 1 _ + J ]3/2 + lw 

which is always positive for any z e [—1, +1] and for any rj > 0. 

Thus, we have only to look for the non zero solutions z of equation 

/-M=0. (53) 

To this end, we consider the function rj(z), defined by ()33[) . which satisfies the 
implicit equation 

f.[z,r}(z)}=0, V«€ (0,1). 

Thus, the inverse function z — z(rj) of rj(z) gives the solutions of equation (|53|) : in 
order to count the branches of the inverse function z = 2(77) we compute the first 
derivative 



rf{z) 



yx+i g{z)- g{-z) 



where 



Since 



[1-Z 2 ]3/2[(1 + Z )<x-(1- Z )<r]2 1 

g(z) = {az 2 - az + + Z y . 



lim 77' (z) = 

then a bifurcation of the stationary solution occurs at z = for 

rf = lim 77(2) = 2 CT /cr . 

Furthermore, a straightforward calculation gives also that 
lim,"(z) = -f(^-3 ff -l) 

z^0+ 3(7 



and 



> if a < a t h 



res 



huhl 



lim 77" (2) < =0 if a = cr threshold (54) 

< if a > Uthresh old 

where 

_ 3 + 713 

^threshold — ^ • 

Hence, we can conclude that in the case a < (^threshold then we have a supercritical 
pitchfork bifurcation at z = (see Fig. Q]- panel (a)), and for a > crthreshold then 
we have a subcritical pitchfork bifurcation at z = (see Fig. [T]- panel (b)). 

Finally, we only have to count the number of branches of the function z(rf) 
and thus we look for the number N of the solutions (counting multiplicity) of the 
equation 

h(z) = 0, h{z)=g{z)-g{-z) , (55) 
for z in the interval z G (— 1, +1). 
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Lemma 4. Let N be the number of solutions z of the equation h(z) = in the 
interval [— 1,+1], counting multiplicity. It follows that z — is a solution with 
multiplicity 3 if a ^ o~threshoid, and with multiplicity 5 if a — a threshold- Further- 
more, it also follows that 

N= { I l f/< threshold (56) 
[5 if a > O 'threshold 

Proof. We may remark that if z* is such that h(z*) = then h(—z*) = 0, too; 
furthermore h(±l) = ±2 CT ^ 0. First of all we see that z ~ is a solution of ([55)1 
with multiplicity 3 for any a ^ crthreshoid] indeed, a straightforward calculation 
gives that 

h(0) = ti(0) = h"(0) = and ti"{0) = Aa(-a 2 + 3ct + 1) . 

Then /i"'(0) ^ if a ^ a threshold- If cr = a threshold then a straightforward 
calculation gives that h"'(0) = h IV (0) = and 

h V (0) = — 86((Tj fcresftoW — ^^threshold + ^ a threshold ~ ^threshold — 6) 

= 24(3 + \/l3)(4 + \/l3) > 0. 
Hence, it follows that 

if CT > (J threshold 

N is { = 5 or > 9 if a = ^threshold 

if a < a thresh old 

where N is number of solutions, counting multiplicity, of equation f—(z,rj) = 
Indeed, we see that 

lim T)(z) = +00 . 
z— >±i 

Then , in the case a > Othreshoid since lrnij^o ff'(z) < then there exists two non- 
zero solutions of equation (|55j) in the interval (—1, +1) at least; hence, the number 
N of solutions of equation ([55]), counting multiplicity, is N > 5. 

In the opposite case a < o 'threshold it follows lim 2 _>.o rf'(z) > 0, then we have two 
cases: or equation ([5"5"j) docs not have solutions z E (—1, +1), z ^ 0, and in this case 
N = 3; or equation (|55|) . counting multiplicity, has other solutions z G (— 1,+1), 
z 7^ 0, and in this case the number of such a solutions is bigger than 4, in this case 
N > 7. 

Finally, in the case a — a threshold it follows that lim^o rj'{z) = lining v" '( z ) 
and 

TV , N 6 • 2 at ^ h ° ld (829^ + 2989) 

lim V IV z = 7 i '- > 0, 

*-*> w 5 (649 + 180VT3) 

hence iV = 5 or N > 9. 

If we can prove that N < 5 then the Theorem is completely proved. 
To this end we set 

1 — z 

V = -r— > V e (°) +°° ■ 
1 + 2 

Hence, equation = in the interval (—1, 1) reduces to the equation of the form 
p a {y) — where 

Pa(y) = y a (y 2 + by + a)-(ay 2 + by + l) 

= y a+2 + by a+1 + ay° - ay 2 - by - 1 
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and where 

a = 1 + 2a, b = 2 - 2cr . 

We remark that if y* > 1 is a root of the polynomial p a corresponding to a give 
z* > 0, then 1/y* < 1 corresponds to — z* and it is a root, too. We remark also 
that: 

- for any a 7^ | [3 + \/13] then p CT has solution j/ = 1 with multiplicity 3; 

- for a = 4 [3 + then p CT has solution y = 1 with multiplicity 5. 
We assume, for a moment, that a is a positive integer and we see that: 

- a = 1: in such a case p\ — (y — l) 3 which has only the solution y = 1 in 
the interval (0, +00) with multiplicity 3; 

- er = 2: in such a case j>2 = (2/ — l) 3 (y + 1) which has only the solution y = 1 
in the interval (0, +00) with multiplicity 3; 

- cr = 3: in such a case P3 — (y — l) 3 (j/ 2 — 2/ + 1) which has only the solution 
y = 1 in the interval (0, +00) with multiplicity 3; 

Now, we are looking for the number N of real solutions, counting multiplicity, 
of the polynomial p a in the interval (0, +00). We already know that for a = 1, 2, 3 
then iV = 3; we also already know that for a = 4, 5, 6, . . . then N > 5 and, in order 
to get an upper estimate of N, we make use of the Budan- Fourier theorem [29 . 

If we denote by v(y) the number of sign changes in the sequence 

then the Budan- Fourier theorem applied to the polynomial p a with degree + 2, 
where cr = 4, 5, 6, . . ., states that 

N < |u(+oo) - u(0)|. 

Since 

lim p^ > 

y— y+00 

for any n = 0, 1, 2, . . . , a then w(+oo) = 0. On the other hand we observe that 



Pa(0) 


= -1< 


pUo) 


= -6 > 


p"(o) 


= -2a < 


p£*>(0) 


= if2<n<cr 


^ CT) (o) 


= ala > 


^ CT+1) (o) 


= 6(ct + 1)!<0 




= a(a + 2)!>0 



since a = 1 + 2<r > and b = 2 — 2cr < for a = 4, 5, 6, . . .. Then, y(0) = 5 and so 
we can conclude that 

N < \v(+oo) - u(0)| = 5 

Therefore, the number of solutions y € (0, +00), y ^ +1, is exactly equal to 2. 

We prove now that N < 5 even for any positive not integer a. In order to prove 
that N < 5 we make use of an extended version of the Budan-Fourier theorem [TU] 
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applied to the polynomial Pa(y) for y G (0, +00). If we assume, for a moment, that 
a > 2 (the case 1 < a < 2 can be similarly treated) then we set 

r = 0, ri = 1, r 2 = 2, r 3 = a, r 4 = a + 1, r 5 = cr + 2 

where (mimicking Example 2 in [10j ) 

a = -1, ai = (2cr - 2), a 2 = -(1 + 2a), a 3 = 1 + 2cr, a 4 = -(2a - 2), a 5 = 1. 
In this case we have 6 functions gj(y), j = 0, 1, 2, 3, 4, 5 such that 

95 (y) =p<j{y) 

and 

lim 9j(y) = +°° 

and 

<7 (0) = 2(<r + l)(<7 + 2) > 
5 i (0) = -2(cr- 1) 2 (<t + 1)ct <0 

52 (y) = (l + 2<7M^l)(a-l)r 3 + O(r 2 )>0a S H0 + 
53 (0) = -2(1 + 2a) < 
54 (0) = 2(ct - 1) > 
ffis(O) = -K0 

Then the sequence gj has zero sign changes at +00 and it has 5 sign changes 
at + , i.e.: v(+oo) = and v(0) — 5. Therefore, Theorem 1 [TU] implies that 
N < \v(+oo) - v(0)\ = 5. 

It remains to consider the case < a < 1. In order to look for the solutions 
y > of equation p a (y) = we observe that these solutions are such that 

a _ ay 2 + by + l 
y 2 + by + a 2 

where the l.h.s. of this equation is a monotone increasing function, while the r.h.s. 
is a monotone decreasing function for < a < 1. Hence, the number of solutions, 
counting multiplicity, of the equation p a (y) = is N = 3. □ 

The proof of the theorem is so completed. 

Remark 12. From Lemma 4 it turns out that when a < a threshold then equation 
r)'(z) — has only solution z — and therefore, under such condition on a, we 
only observe a bifurcation of the stationary solution at \n\ = rj* . On the other side, 
when a > athreshoid then the number of solutions ( counting multiplicity) of equation 
r)'(z) — is 5, since the solution z = has multiplicity 3 then the other 2 solutions 
are ±z + , where z + £ (0,1), and they are associated to saddle points appearing at 
\n\ = r/ + , where rj + = r)(z + ). 

Remark 13. We just point out that in the case ofrj < then we can apply the same 
arguments; we only have to emphasize that for negative values of r\ then equation 
f-(z,rj) = does not have non zero solutions and that bifurcations come from 
equation f + (z,rj) = 0. 
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Figure 2. In this figure we plot the graph of the values of the 
function E versus the nonlinearity parameter rj for nonlinearity a = 
1 < ^threshold- For 7] = ±77*, i]* = 2 for a = 1 , a bifurcation occurs 
and a new branch corresponding to the asymmetrical stationary 
state appears. Line (s) denotes the symmetric stationary solutions, 
line (a) denotes the antisymmetric stationary solutions, and (as) 
denote the asymmetrical stationary solutions. 



Remark 14. For large a the roots y < 1 of the polynomial p a (y) are asymptotically 
given by the roots of equation 

(1 + 2a)y 2 + (2 - 2a)y + 1 = 0. 

That is 

Hence, the solution z + of equation rj ' (z) = is asymptotically given by 

a cr z 

and we have that 

r/+ = V2ea [l + 0(cr~ 1 )] 

in the limit of large a. 

Remark 15. The frequency A of stationary solutions of equation 125\) are thus 
given by 

X = Q, + ojE 

where E = E(z) is the multivalued function given by i51\) , where z — z(rf) are the 
roots of the equation f±(z) = 0. For the graph of the functions E{z), depending on 
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FiGURE 3. In this figure we plot the graph of the values of E as 
function of the nonlinearity parameter rj for critical nonlinearity 

d = ^threshold- 

T), we refer to the Fig. [H Fig. and Fig. [7J We observe the following behaviors 
(where we assume 77 < for argument's sake): 

- When -rf < r/ < for a < a thresho i d , or -77+ < r\ < for a > ^threshold, 
then we only have the linear stationary states. 

- When rj < —rf and a < crthreshoid, then the symmetric solution bifurcates 
at rj = —rj* and then we have 4 stationary solutions: the two linear sta- 
tionary states and two new asymmetrical stationary states; a similar picture 
actually occurs also when a > crthreshoid, but in this case the two new asym- 
metrical stationary solutions don't come by a bifurcation of the symmetric 
stationary solution, but they come from a branch of saddle points. 

- When ~rf < rj < —77+ and a > crthreshoid, then a couple of saddle points 
occurs and thus we have 4 asymmetrical stationary solutions. Two of them, 
denoted as ( asl ), are much more localized on a single well than the ones 
denoted by (as2). 



4. Dynamical stability 

The time-dependent equation (|24p , when projected on the one-dimensional spaces 
spanned by the single- well states ipn and tfz, and on the space LT c L 2 (R d ), takes the 
form 

ia' R = -a L + r R 

ia' L = -or + r L (57) 
W c = L[H -n]Tp c + r c 
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Figure 4. In this figure we plot the graph of the values of the 
function E versus the nonlinearity parameter 7/ for nonlinearity 
(7 = 5 > (^threshold- At |?7| = 7? + , 7/+ ss 4.41 for a = 5, a couple 
of saddle nodes appear, and the corresponding branches, denoted 
(asl) and (as2), are associated to asymmetrical stationary solu- 
tions; asymmetrical solution (as2) then disappears at \q\ = 77*, 
rf = 6.4 for a = 5. 



where we have set tp — > e^ i ^ T ' 0J ijj(x, r). We call two-level approximation the 
system of differential equations coming from (|57p taking ip c — and neglecting the 
exponential remainder term in r R ^(a R , cil, 0) (see Lemma [5]); in such a case the 
two-level approximation takes the form 

ia' R = -a L + i]\a R \ 2a a R , |2 , |2 , 

,•„/ _ „ , ' „ \2a n , \ a R\ + ml =1 (58 
%a L — — a R + r]\aL\ ol 



We may remark that the two-level system (|58|) takes the Hamiltonian form 

iA' = d A H, A = (a R ,a L ) , 
with Hamiltonian function 



H 



(a R a L + a L a R ) - JL- (\a R \^ + \a L \^ +1 A 
a + 1 V / 



(59) 



corresponding to the energy functional restricted to the two-dimensional space 
spanned by the two single-well states. The stationary solutions of the two-level 
system ((58)) are associated to stationary points of the energy functional W, then we 
can attribute them some stability /instability properties in the sense of the theory 
of dynamical system. In particular, let 9 — &rg(a R ) — arg(az,) be the difference be- 
tween the phases of a R and cil, and let z = \a R \ 2 — |ol| be the imbalance function, 
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then system (|58|) takes the Hamiltonian form 

e = d z n 

z = —deH 
where the Hamiltonian (|59p takes now the form 



(60) 



z A cos ( 



n 



i 



1 i \ <j + 1 / 1 \ cr + 1 

1 + Z \ I 1 — z 



In order to study the stability properties of the stationary solutions of equation 
(f6"0"]) we have to consider the matrix 



Hess = 



a-u 

ae 2 




at the stationary points. Since the trace of Hess is zero then we have that the 
stationary point is a circle if det Hess > 0, and it is a saddle point if det Hess < 0. 

4.1. Dynamical stability of the symmetric and antisymmetric stationary 
states. We consider, at first, the symmetric and antisymmetric stationary states 
corresponding to 8 = and z — (symmetric), and 9 = ir and z = (antisymmet- 
ric). A straightforward calculation gives that 

det Hess\ e=Q z=0 = 1 + 7]^ and det Hess\ g=7r ^ Z=Q = 1 - 77^ . 

Then, it follows that the symmetric stationary solution is dynamically stable for 
any 77 > —77*, and the antisymmetric stationary solution is dynamically stable for 
any r\ < if , where 77* = 2 <T /cr. 

4.2. Dynamical stability of the asymmetrical stationary solutions. For 

argument's sake let us assume 77 < 0. Then the symmetric stationary solution 
bifurcates and new asymmetrical solutions appear, they correspond to 6 = and 
the values of z are the non zero solutions of the equation f + {z,rj) — (in fact, 
we have assumed 77 < 0; in the case of 77 > 0, as considered in $3] for the sake of 
dcfiniteness, then the stationary solutions corresponds to the roots z of equation 
f-(z,rj) = 0). A straightforward calculation gives that 



det Hess\ g _ Q — 



(l-z 2 )- 




(7—1 / ^ \ (T — 1 



By the relation 77 = 77(2) implicitly defined by the equation f+(z,r)) — it follows 
that 

det Hessl - 9{z) - g( ~ z) 

le =°< "="( 2 ) (1 - z 2 ) [(1 + z)° - (1 - zY] 

where it has been already proved that the equation g(z) — g(—z) = has a solution 
at z = with multiplicity 3 (multiplicity 5 if a — o 'threshold)- Since this equation 
has no other solution for a < athreshoid, since q(z) = q(—z) and since 

lim_ det Hess\ g= ^ v ^ [z) = +00 

then 

det Hess\ e=0> n=v(z} > , Vz ^ 0. 
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Then, the asymmetrical solutions, if there, are stable. On the other side, for 
a > athreshoid then the equation g(z) — g(—z) — has three distinct solutions; 
hence, by means of the same arguments as before, it follows that the branch (as2) 
is dynamically unstable and the branch (asl) is dynamically stable. 
We can collect all these results as follows (see also Fig. []}. 

Theorem 2. Let us consider the stationary solutions of the two level approxi- 
mation \58\) that coincide, up to an exponentially small term, with the solutions 
given in Theorem]]^ The symmetric and antisymmetric solutions of the two-level 
approximation are such that: 

- for any a > 0, the symmetric stationary solution (s) is stable for any 
t] > —if , and it is unstable for any rj < —rf; 

- for any a > 0, the antisymmetric stationary solution (a) is stable for any 
77 < 77* 7 and it is unstable for any r] > rj* . 

The asymmetrical solutions of the two-level approximation are such that: 

- for any a < athreshoid the asymmetrical stationary solution (as) is stable; 

- for any a > a t h. r eshoid the branch (as2) of the asymmetrical stationary 
solution there exists for any n + < \n\ < 77* and it is unstable, the other 
branch (asl) of the asymmetrical stationary solution there exists for any 
n + < \r]\ and it is stable. 



5. Orbital stability 

In this section our aim is to study the orbital stability of the stationary solutions 
of the NLS (p}. So far we have considered both cases of attractive and repulsive 
nonlinearity for any couple of eigenvalues A± . Hereafter we consider only the first 
two eigenvalues and we assume to be in the attractive nonlinearity, that is: 

Hypothesis 4. LetX± be the first two eigenvalues of Ho- Letn= ^(ip' T ^ rl , gip^ 1 ) 
be the effective nonlinearity parameter in Eq. \lJ$ where (ip 1 ^ 1 , gip^ 1 ) > 0; we as- 
sume that 

e < that is 77 < . 

If we rescale the solution if) as <f> = \e\ 1 / 2cr ip, then equation (jTJ) is equivalent to 
the equation 

i^=#o0-<#| 2< ^, U\\ = N 1/2CT - (61) 
The stationary solutions of the equation 

Ho4>x,e-9\^A 2a ^e-X(t>^ =0, X = n + ujE, (62) 

are associated, by means of the scaling, to the stationary solutions ip s E , if>% and 
ip'jz given in Theorem Q] where E = E(e) is a multivalued function and where the 
stationary solutions are now denoted by 



A,e 
a 

X,e 
as 
X,e 



iasl and J - as2 



symmetric stationary solution 
antisymmetric stationary solution 
asymmetrical stationary solution for a < ^threshold 
asymmetrical stationary solutions for a > athreshoid 
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If we consider a general stationary state, we denote the solution by (f>x,e and ipE, but 
if we want to distinguish the branches, we insist, in such above way, by denoting s, 
a, as, asl and as2, on each shoulder of solutions. 

Here, we consider the orbital stability for the symmetric stationary solution 
4>\ e and for the asymmetrical stationary solutions <p^ s e that bifurcate from the 
symmetric one. 

Definition 1. The family of nonlinear bound states {e la 4>\. e , a € K} is said to be 
orbitally stable in H 1 (M. d ) if for any k > there exists a S > such that if (j>o 
satisfies 

inf Uq - e ia cl> x , e \\ H i < 8, (63) 

then for all t > 0, the solution cj)(t) of H61\) with 0(0) = (j>o exists and satisfies 
inf ||0(-,;)-e la 0A,e||ffi <K. 

Otherwise, it is said to be unstable in if^R^). 

The main result of this section is the following: 

Theorem 3. Fix any H > be sufficiently small such that H £ (0, /13) for some 
/13 > small enough. Then, the following statements hold. 

• Let a < cr t hreshoid- The symmetric solution corresponding to z s = 0(e~ p / h ) 
is orbitally stable in H 1 for \ri\ < rf . At the bifurcation point ?/ = rf , there 
is an exchange of stability, that is, for > rf , the asymmetric solution is 
stable in H 1 and the symmetric solution is unstable. 

• Let a > a threshold- By Theorem\]\ two couples of new asymmetric station- 
ary states, denoted by ip asl and ip as2 appears at \n\ = r; + . For \rj\ > r] + , 
ip asl is orbitally stable in H , ip as2 is unstable. On the other hand, the 
symmetric state is orbitally stable in H 1 for \n\ < rf , and unstable for 
\V\ > rf. 

As a standard method to prove the orbital stability of a stationary solution <f>\^, 
the following proposition is well known. We first define L^ e and L^ e , which are 
respectively the real and the imaginary part of the linearized operators around a 
real valued stationary solution 4>\ i(i : 

L x / = L+[<j> Xie ] = H - A - (2a + l)g\<f>x,e\ 2a , 

L x _: e = L^ c ]=H -\-g\ct>xA 2a - 

It is clear that L^' e 0A,e = since 4>\ t c is a solution of (|6"!2)) . Moreover, L ,' e and L x _f 
are self-adjoint operators on L 2 (M. d ) with domain H 2 (R d ). The essential spectrum 
of these two operators coincides with the interval [V^ — A, 00) with — X > 0, 
since 4>\ >t vanishes at infinity; indeed, V is bounded, and we can apply the proof 
of Theorem 1 in [15] . regarding the term V<f>\ ye of (|62]> as one of nonlinear parts. 
There are also finitely many of discrete spectrum and ad(L^ e ) C (— 00, — A) 
(see 01). 

In order to prove the orbital stability we make use of the following criteria (see, 
e.g., [H] or Part I of [19]). 

Proposition 1. Suppose that L^ is nonnegative. Let F(X) — ||0a. c || 2 . 
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(1) If L+ e has only one negative eigenvalue, anddF/dX < 0, then, <f>\^ is stable 
in H\R d ). 

(2) If L^ e has only one negative eigenvalue, and dF/dX > 0, then, (f>\^ is 
unstable in H l {R d ). 

(3) If L+ e has at least two negative eigenvalues, then, <p\ te is unstable in 

Remark 16. For the instability (3), it is enough to find a vector p G H 1 such that 

(L$ e p,p) <0, pl^mi 2 . (64) 

(see for e.g., [71 118] ). As we will see below, "L^ e is nonnegative and L^! e has two 
negative eigwnvalues" occurs only for the symmetric stationary solution <p\ t . In 

this case, we can find the normalized antisymmetric solution \\ as the vector p 
satisfying the property {6$ for h small. 

We shall therefore check the following properties: 

• the number of negative eigenvalues of L+ e ; 

• L^' e is a nonnegative operator; 

• (Slope condition) the sign of the function dF(X)/dX. 

5.1. Number of negative eigenvalues of L^ e . First we consider the number of 
negative eigenvalues of L+ e . We will prove that: 

Lemma 5. Let h* > small enough as in Theorem^ there exists hi G (0,fi*) 
such that for any h G (0, hi) the following statements are satisfied. 

(i) Let X be the energy level associated to the symmetric stationary state <f>\^ — 
4>\ e . Then, L^ e admits only one negative eigenvalue provided that |^| < 

rf . On the other hand, L^! e admits two negative eigenvalues provided that 
\n\ > rf. 

(ii) Let X be the energy level associated to the asymmetrical stationary state 

(f>\,e = (t>\% if a < ^threshold, and 0A,e = 4>T,e and 4>\,e = <f>T,e > */ a > 

^threshold- Then, L^ admits only one negative eigenvalue. 
Proof. We set 

0a, £ = o£Vfl + o^'Vl + 4>Y, \4T + I«l'T + UYf = 

H^'ll = lel^l^dl, ^c = *i-(ak + 4)i 

where tf>E is a stationary solution obtained in Theorem [TJ 

We consider the eigenvalue problem L+ e u — with u G H 2 (R d ) and where 

\{iuj\ < Ch 2 . (65) 

By setting u = dR(p R + CLLfL + u c with u c G Tl c L 2 , then the eigenvalue problem 
takes the following form 

u>Lia R = a R Vt - a L uj - Xa R - (2cr + l)((p Rl g\(j>\^\ 2 ' y u) 
LU[ia L = a L tt- a R uj - Xa L - (2a + l)(<p L ,g\(j)\^\ 2 ' y u) . (66) 
uj[iu c = (H a - X)u c -II c (2o- + l)g\4> x ,e\ 2 ' 7 u 
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The last equation reads as 

I — [Ho - A - c^p 1 n c (2a + l) 5 |0A, £ r] «c 

= (flo - A - cj^) _1 n c (2cr + l)y|(/!)A,e| 2cT (ai?^ J R + ol^l) 
Since Ho — A > Cft, when restricted to H C L 2 , and since (1551) . then we have 

|| (h - A - w^)^ 1 n c || £(i 2^ ff 2) < dft -1 . 

Here, we recall that, from ([20)) and (|36]l . 



with ao = 1 
small ft, 



d(2er-l) 
4 



| 2CT || = k|||9|V B (e)| 2<7 |l<CN^" a ° 

Thus, if ojfi = 0(ft 2 ), we get from (fl7|) . for sufficiently 



|| (flo - A - c^/i)- 1 n c (2a + l)<#A,eHU(£^ 2 ) < C 2 (2<t + l)|e|ft- Q » < 
Namely, if /i satisfies the condition ([r?5|) then the inverse of the operator 

J - [flo - A - ujfi]- 1 n c (2a + 1) 5 |^a, £ | 2<t 
exists. Accordingly, the third equation in (|66j) has a solution 
m c := u c (fi,X) 

where 

|2(T 



Q[/x,< 



J-(flb-A-wM) 1 n c (2(T + l) 5 |^ 



x (H - A - w/i) 1 n c (2a + l)<# A ,e| 2CT : i 2 ( 



*) -> H 2 (R d ), 



and 



IIQ[^0A, £ ]||£(L^) <a|e|^ Q0 . 
The bound C a is uniform in ft on D, where D is defined in Lemma [31 and for any 
fi such that \fMu\ < Ch 2 . In fact, by the same arguments the same estimate holds 
true also for the derivative of Q with respect to /i: 

dQ 



dfi 



< C|e|ft- Q « (67) 

C(L 2 ^H 2 ) 

for some a' Q > 0. We insert this expression of u c into the system (|66p. and we have 
Lyapunov-Schmidt reduction of (j6"6")l as follows. 

ujfiCLR — a R Q — a^uj — Xa R 

-(2d + l)(<p R ,g\<f> Xte \ 2a (I + Q(ji, 4>x, e ))(a R <p H + a L cp L )), 
u)/j,a,L = a,i,Q, — a R u) — Xcll 

-(2(7 + l)(<^, 5 |<^ e | 2CT (7 + Q(p, (j>x, e ))(a R <p R + a L <p L )). 

This system can be rewritten under the following form. 

,0, 



(N + fil- vC) 



a R 
ah 

where we recall that A = f2 + ojE and where 

'1 



(68) 



N 



a, 1 
1, P 



I 



1 



c 



Ci, C*2 

C3) C4 



a = E + (2a + 1)M|4| 2(T , /3 = J5 + (2(7 + l)M|a£ 



A |2o 
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and v = e lp l h for any p' <E (0, p) as in (|45j) . For v ^ 0, we have put 
d = d(4, al E, n;h) = - {2U + 1)|e| { (<p R) 9{\M 2 ° - I4y fl | 2 >fl) 

+(^,,3|V'e| 2c7 Q(^,0a, £ )^)}, 

C 2 (a x , a x ,E, W ^ = - (2<T + 1)|e| {(^, g|^| 2CT (/ + Q(a*, <M)^>}> 



C 2 

C*3 = C2, 



a 



4 



C4(o^ai^ w n) = - (2CT + 1)|e| {(^ ;g (|v^r - I«V| 2ct Vl) 

+ (<PL,g\4>E\ 2a Q(V,^\,e)VL)}- 

If ^ = 0, then p, are the eigenvalues of N and they are the solutions of the 
equation 

PQi) = fJ, 2 + (a + (3)p + a/3 - 1 = 0, 

which always has only two real different solutions /ii, [ii since (a + /3) 2 — 4a/3 + 4 = 
(a — /3) 2 + 4 > 0. In particular, these two real eigenvalues are both negative or 
both positive if a/3 > 1, or only one is negative in a/3 < 1. 

To investigate the sign of a/3 — 1, we consider, at first, the case of the symmetric 
stationary solution corresponding to z x = z s — (see Theorem Q] and Remark 
[7|). Then (hereafter, for the sake of simplicity, we denote by ~ that we have an 
exponentially small term) 

4 = a£ = i=, £~-l-|ry|i 

and 

a = P = E + \ V \^±±~-l + lr ] \^ 

Hence, condition a/3 > 1 is equivalent to the condition \r)\ > 77* = — (and in 
such a case both solutions are negative), and condition a/3 < 1 is equivalent to the 
condition \rj\ < rj* = — ; provided h is small enough. 

We consider next the case of the asymmetrical stationary solution corresponding 
to z x 7^ 0. In such a case we set a — \\rj\(p x ) 2<J and b = \\'q\(q x ) 2 ' 7 , then 

a ~ E + 2(2cr + l)o , P~E+ 2{2a + \)b 

and 

E ~ -ijl-(z^- 2[a(p X ) 2 + b(q X f] 
Hence, condition a/3 < 1 is equivalent to the condition 

l{z x ,a) -1<0 

where 



£(z,a) := y/l-z 2 + 2 [(p 2 -2a- I) a + q 2 b] 



Vl-z 2 + 2 [ap 2 + b(q 2 -2a- 1)] 
[(-1 + z + Aza)(l + zf + (1 - z) a+ ^ 



(1- z 2 )[(l- zY -(1 + z)-] 2 



[(1 + z + iza)(l - zf - (1 + z) 
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since 

P ~ \J~Y~ and q r 
Then, a straightforward calculation gives that 

(l- z 2 )[{i + z y z y\ 2 



-•- 1 -<*- 1 >{^ +!, * +j '4t$ 

4zcr(l + z) 2ct 



(l-z 2 )[(l + z) CT -(l-z) CT ] 2 



X - A_ [-1 - y 2ff+1 + (1 + 2a)y CT - (1 + 2 ( 7)y CT + 1 ] 

1 — z 

where we have set y = € [0,1]- We then consider the sign of the following 

1 ~\~ z 

polynomial in the right hand side above. 

q(y) := -1 + y 2a+1 + (1 + 2a)y° - (1 + 2a)y CT+1 . 

It is in fact easy to conclude that q(y) < for any y € [0, 1]. Indeed, 

q(y) < -1 + (1 + 2a)y a - (1 + 2a)y a+1 < 0. 

Now, we wish to investigate the sign of eigenvalues for the case v ^ 0. Recall 
that the effective nonlinearity parameter r\ satisfies \rj\ < C for some constant 
C > 0. Also there exist ho € (0, h*), and a compact interval i^ft such that the 
two eigenvalues of the matrix N, 

Vi = \{-(a + P)- v/(a-/3) 2 +4}, M2 = i{ -(a + /3) + ^(a - P) 2 + 4} 

belong to K% for any ft £ (0, fro)- Then we see that Cj = Cj(a^,a^,E,fi,h) 
are bounded, together with their first derivatives, on D x K% uniformly for any 
H £ (0, Hq): indeed, there exists a constant C > such that 



and this right hand side is bounded because of (fl4"|) . (f2"0)l and It also follows 

that if 1 < 2cr, 

^- 1 |(^, 5 (|^| 2 ' 7 - |aW| 2CT )^)l < v^CQIvrVlW + ||Vc||)(l + 

whose right hand side is bounded, noting (|20] ) . (f2"T]l . (|35] ) and (|33 | . If < 2a < 1, 

|(^, ff (|^| 2CT -|a^| 2CT )^)|<C J i&lgWalipL+licfdx 

<C\\g\\ L ~ J ifl^dx + \\9\\L~\a x L \ 2 ° J <f? R <f%dx. 
The first integral is estimated as follows 

/ vim 2 °dx< \wi\\ L , ^\m 2a \\L^m\ 2a \w R \\iv^ m 
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by means of the Holder inequality, where q = i > 2 and p — jzrz- Inequalities 
([Ml) and (fH)) yield that this right hand side is exponentially small. Similarly, the 
estimate of the second integral follows 



for some a > 0. As for the derivatives of Cj, the analyticity in /i of (Hq — A — w/i) -1 
ensures their regularity, and the uniform boundness follows from ()67[) . 

We come back to the problem (|68|). This problem is mapped to the problem to 
find the roots of the following characteristic equation. 

D(a R: a L ,E, fi, v) = det(N + \il - C) = 0. 

Concretely, 

det(iV + /xI-C) = (a + M-^Ci)( i 8 + j u-^C4)-(l-^C 2 )(l-i/C , 3) 

= fj 2 + {(a + j8) - i/(£7i + C 4 )}^ + a/3 - 1 

-f (C 2 + C 3 + aC 4 + + i/ 2 (Ci<7 4 - C 2 C 3 ). 

Putting 5(m, i>) = -(Ca + C 4 )m - (C 2 + C* 3 + aC 4 + /3C*i) + v(C x C A - C* 2 C 3 ), we 
have 

D(fi, v) = P{n) - vS{n, v) = 0. 

We note that by the above arguments, 5(/x, v) and d^P{^) is uniformly bounded 
on D x if^,, for any h £ (0,ho). It is also seen that D(fi,v) is a C 1 function in 

D(/ii,0)=D(/ia,0)=0, 

By applying Implicit Function Theorem, there exist £o > such that there exist 
two real solutions and of £>(jU, = for \v\ < £o and that 

WM = /ii-«4^ + 0(^), (70) 

i*iy) - ^ 2 -4^ + <V). (7i) 

(72) 

Therefore, for any e > there exists h\ £ (0, So) such that \n\{v) — H\\ < e, and 
that |/x 2 (f) — /i 2 | < e and /ii(^), jU 2 (f) € -fOi f° r an Y ^ G (0, fi-i). We remark 
here that L + ' e has at least one negative eigenvalue since (L + ' e 0e,A, 0c,a) < 0. As 
a consequence, for the symmetric solutions, L^' e has two negative eigenvalues if 
1 77 1 > rf and has only one negative eigenvalue if |ry| < rf. For the asymmetric 
solution, _L + ' e has only one negative eigenvalue. The proof of Lemma [5] has been 
completed. □ 

5.2. i^' e is a non-negative operator. Next our aim is proving that L x _f has 
no negative eigenvalues. Since the symmetric solution ipg, i.e. 4> s e \ 1 is positive by 
means of a suitable choice of the phase, L^ e [(f> s e ,] is non-negative. However, we 
do not know the sign of the asymmetric solutions and we repeat here the same 
argument as in Lemma [5] for L^ e . 
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Lemma 6. Let (f>x,e be the symmetric and asymmetrical stationary solution asso- 
ciated to the level X. Then there exists H2 € (0, Hi), where Hi has been defined 
in Lemma\5[ such that for any h 6 (0, hi), Lj has no negative eigenvalues, more 
precisely, L x _f has a zero eigenvalue and one positive eigenvalue uifi = 0(H 2 ). 

Proof. The eigenvalue problem L^ e u — (w/i)tt with u £ H 2 (M. d ), where \ojfi\ < Cti 2 , 
takes the form 

iop,a R = a R n - a L u - Xa R - a R (tp R , gl^x.^u) 
u>p,a L = a L Vl- a R uj - \a L - a L (Lp L ,g\4>\. t \ 2(I u) 
uj[lu c = (H - X)u c - IlcglfixJ^u 

where we put u — a R tp R + ol^Pl + u c , u c € II C L 2 . We remind that /i = 
is a solution of the eigenvalue problem since L^tfix^ — 0, we then apply again 
the same Lyapunov-Schmidt reduction as in Lemma [5] in order to compute the 
other eigenvalues of L^ e such that \/j,oj\ < CH 2 . This eigenvalue problem can be 
rewritten, assuming \cufi\ < CH 2 , as follows. 



(N' vC) 



a R \ (0 



\-i _ i a ' i M r" — f^'ii C'2 1 fi _ /-•/ 

JV — 1 qi I ' u — \ /~<i nl ) °3 — °2i 



where 

N' = ( , , - , 

a' = is + H|a^r, /3' = £+M| a £| 2CT . 

Remind that is defined in Lemma As in the proof of Lemma O it suffices 
to know the sign of a' ' fi' — 1 . We compute the case of the asymmetric solutions 
corresponding to z A 7^ (in the case of the symmetric solution corresponding to 
z x = we follow the same arguments). In this case, 

ol v /l-(z^-| 7? |{(/) 2 -+ 2 + (q X ) 2 °+ 2 } + \ V \(p X r, 

fi' x /l (z x ) 2 \ V \{(p X ) 2a+2 + {q X ? a+2 } + M(q X ) 2 \ 

-1 2z x 

x 



r/l + z A \ <T (\-z > 



By direct computations it is not difficult to obtain that 

z x - 1 , z x + l 



Therefore, a'f3' ~ 1 and 

a' + fi' 



which implies /ii/i2 ~ and [ix + > for the eigenvalues of iV 7 . We may 
assume without generality that |/ii| is very small and \x-i is positive. It follows 
from the same arguments as in Lemma [5] that the perturbed matrix N' — vC has 
two different eigenvalues /ii(^) and ^(v) verifying (|70|) and (|7ip. Since we know 
that Lj e has always zero eigenvalue, and perturbed eigenvalues are continuous with 
respect to v, we conclude that \i\{y) = and \ii(y) > 0. □ 
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5.3. Slope condition. In order to check the slope condition, we consider the fol- 
lowing quantity. 

F(X) = ||0 A , e || 2 = 



and we remark that 

dF(X) 
dX 



dX 
~dl 



l/a 



\(l-a)/a 



1 



C R a' 



dc 



dE 
drj 



-l 



dE 



dc 



Thus, we only have to check the sign of ^ for the symmetric and asymmetrical 
stationary solutions. 



5.3.1. Estimate of the stationary solutions as function of the non-linearity param- 
eter. The stationary solution 

ijj = a R ip R + a L tp L + ip c , 

of equation (P?5l) associated to the energy level E depends on the value of the 
nonlinearity parameter r\ — ec/u, where c = C R — C'l is defined in equation (|30|) . 

In particular, in Theorem 1 we have proved that, locally, there is a correspon- 
dence one-to-one from 77 to the solution p, q, a, /3 and E (up to the gauge choice 
of the phase, where we set 8 = a — (3) of equation (|47|) and tp c of equation (f37|) : 
provided that rj ^ ±n + and rj ^ ±77*. 

In order to see the sign of dE/drj, we wish to obtain the estimate of the first 
derivative of p, q, a, ft, E and ip c as function of rj. To this end, let 

D" = {(p, q, a, p, E) € [0, l] 2 x [0, 2tt) 2 x K : p 2 + q 2 < 1, \uE\ < Cfr 2 } 

for some C > fixed; and let 

$ : Rx H 2 x D" -> H 2 x R 4 
(rj,tp c ,p,q,a,l3,E) i-> (F(ip c ),G) 

where F(ip c ) is defined by PU|) and where G is defined by (|47l) with e replaced by 
ur)/c; F(<ip c ) = F(n,ip c ,p,q,a,P,E), and G = G{r},p,q,a,P,E) = (G 1 ,G 2 ,G 3 ,G 4 ). 
For simplicity, we set y = (ip c ,p,q,a, /?, E) G H 2 x D". Since the mapping 



|%,.) :H 2 xD"^H 2 
dy 



is one-to-one at a point n =^ ±77* , ±77+ , we ob- 



tain the unique solution y = 2/(77) of equation = (up to the gauge choice 

of the phase). Furthermore, $(r),y) is C , so the solution 7/(77) is C 1 except for 
rj 7^ ±77*, ±?7 + , and we have 



— + — y' = . 

an ay 



(73) 



Here, ' denotes the derivative with respect to 77, and we use this notation hereafter, 

<9<i> 

too. We will in fact see that -^—(77,-) is one-to-one for any 77 ^ ±?7*,±?7 + in the 
__ y 

proof of Lemma [8] below. Therefore we do not mention the details about this fact 
here. 



The first equation of (|75|) takes the form 
uiEtp' c + uE'tpc = [H 



fi]^ + -n c 7j + ^-U c gW {a' R ip R + a' L <p L + O , (74) 
c c 
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where v = g\ip\ 2a 'ip7 an d 

W=[( C 7 + l)|^| 2<T + ^ 2 |V| 2(ff " 1) 7l, Tu:=u; 

actually, the stationary solution is a real valued function by means of a gauge choice 
(see Remark 7). 

In order to write the other equations of (|73j) we make use of the first two equations 
of (1281) and of the normalization condition: 



Ea R = -a L + ^(ip R ,v) 

l { l PL,v) 



Ea L 

\ a R\ 



-an 

M 2 



(75) 



Now, we get the estimate of the derivative of ip c in Lemma [7] and then the 
estimate of the derivative of p, q, a, j3 and E in Lemma [51 

Lemma 7. Let (a R ,a.L, E) £ D and let r\ satisfying Hyp. let ip c be the solution 
of equation j37\ l. Then 



dip c 



di) 



H 2 



1 + max 



da R 



da L 



dr) 



dE 



dr) 



b(e- p / h ^j as h^O. (76) 



Proof. Since equation (|74[) can be written as 



[(H -n-uE)n c + en c W] i>' c = w£% - eH c gW (a' R <p R + a' L ip L ) + —n cV . 



then 



tp' c = I + (H -Q- wE)~ 



eU n W 



[H - n-ojE}" 



uiE'ipc - eU c gW {a' R ip R + a' L <p L ) + — U c v 



and, by making use of the same ideas applied in the proof of Lemma [SJ it turn out 
that the inverse operator is bounded and ([75]) follows. □ 

Lemma 8. Let \r)\ ^ rj* and \rj\ ^ r) + . Then 



dp 
drj 



dq 
drj 



da 



dy 



dp 



dij 



dE 



dr) 



< C 



(77) 



for some C > 0. 



Proof. Now, in order to give an estimate of the derivative of p, q, a, /3 and E 
we write down the corresponding equations of (|73p . that is we have to consider 
the derivate of equations ([75]). We assume, for the sake of definiteness, that the 
stationary solution corresponds to 9 — (that is ip is a symmetric or asymmetrical 
stationary solution). In fact, we rewrite a R = pe ta and ol = qe l " by means of p, 
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q, a and /3 (where we set = a — f3); so that equation (|75)l takes the form 



' Ep + qcosO- -§^[{tp R ,v)e- ta ] = 
qsm9+^[{ip R ,v)e- ia ] = 

Eq+pcos9~ JLSi [(ip L ,v)e-^] = 




I p 2 + q 2 + 1 1 -0c I 



= 1 



We take now the derivative of both sides with respect to 77, obtaining that 



' E'p + Ep' + q' cos 9 -q6' sin 9 - [(<p R , v')e~ la - ia'((p R , v)e~ ta ] = [(ip R ,v)e-' la ] 

g / sin0 + g0 / cos + ^3 [(<p R , v') e - ia - ia'(<p R , v)e~ la ] = -^3 [{ip R ,v)e~ ia ] 
< E'q + Eq' + p' cos6 - p6' sm6 - ^3? [(^ L , v')e~ i/3 - i/3'{<p L , v)e~^) = ^9? [(tp L , u) e - ,;/3 ] 
-p'sin0-p0'cos0+ [(y> L , u')e"^ - i/3'(^ L , w)e- l/3 ] = -^3 [(ip L , v)e~ lfi } 



We remark that 



{ V R,g\VR\ 2 >R)\a R \ 2 °a R + 0(v) = C RP 2 ° +1 e"* + 0{v) 
(VL,9\VL\ 2 °<PL)\a L \ 2 °a L + 0{v) = C L q 2 ° +1 e^ + 0{u) 
= (<p R ,gW<p R )(p' + pia')e la + q' 0{v) + f3' 0{v) + 0{v) 
= (^,j%)(?' + #')^+p'OM+«'OM + 0{y) 



where 



(<PR, gWifi R ) = C R [(a + 1) + ae aa ] p 2 ° + 0{v) 
(<p L , gW VL ) = C L [(a + 1) + ae^] q 2 ° + ~0{ V ) 
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Therefore, the above system takes the form (where the asymptotics ~ means that 
the remainder term is of order 0{v)) 

( E'p + Ep' + q' cos 9 - qff sin 9+ 

-C^ + !) + ve 2la ]p 2a {p' + ipot) - ia'C R p 2a+1 ] ~ p 2 ° +1 

q' sin 9 + q& cos 9+ 
+ ^3 [C R [(a + 1) + ae 2la ]p 2a {p' + ipot) - ia'C R p 2 ° +1 ] ~ 



E'q + Eq' + p' cos 9 - p6" sin 9+ 

\c R [{a + 1) + ae 2i P}q 2 °(q' + iq$) - ifi'C R q 2 - +1 } 

—p' sin 9 — p9' cos 9+ 
+ zk^ [°Ri( a + !) + V CT ((z' + ig/?') - ;/?'Ch<? 2ct+1 ] 

. 2pp' + 2qq'~0 
that is 



,,2er+l 



M(l + 0(i/)) 



P' 

4 

a' 
V P' J 



where 



M 



( p E- n[{o + 1) + a cos(2a)]p 2 
rj<jp 2<T sin(2o!) 
cos 9 
— sin 6* 
2p 

—q sin 9 + r\ap 2(T+1 sin(2a) 
gcos6» + i]p 2a+1 a[l + cos(2a)] 
— p sin 9 
—p cos 9 






q 
o 

V o 



/ p 2a+l \ 



f T1 




costf 
sin 6» 

£ - ri[(er + 1) + cr cos(2/3)]g 2lT 
r\aq 2a sin(2/3) 
2q 

+q sin 6* 
— qcos9 
+psin9 + r]aq 2a+1 sin(2/3) 
pcos6» + m 2a+1 (j[l + cos(2/3)] 




/ 



We consider now, separately, the symmetric and asymmetrical solutions. 

Symmetric solution. In the case of the symmetric solution where 9 = we 
can choose the common phase a = /3 = 0, by means of a gauge choice. Since 
p = q = ^7=, then a straightforward calculation gives that the matrix M takes the 
form 

det(M) = -8arj2- a (l + 

Hence, for 77 < then det (M) 7^ provided that |7y| 7^ ij* . Hence, we have that 
(1771) holds true. 

Asymmetrical solution. In the case of the asymmetrical solution corresponding 
to 77 < then 9 = 0, we can still choose the common phase a = /3 = by means of 

a gauge choice, and P = \ an< ^ 1 = V ^~T~ sa tisfy equation f+(z, rf) = 0. Then 
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we can set 



2z 



and 



E 



Vi - z 2 + ■ 



l + z 



l + z 



1 - z 



1 - Z 



cr+1 



By means of a straightforward computation it turns out that 



dct M = 8az2 a+1 ^ 



-(4<r + 2)z(l - z 2 Y + (1 + z) 2cr+1 - (1 - z) 2a+1 ][h{z)} 



(1 - z 2 )[(l + z y - (1 - zYY> 

where h(z) — g(z) — g(—z) enters in the definition of rf (see equation (51)). If we 
remark that the function 

Q{z) := [-(4a + 2)z(l - z 2 f + (1 + zf° +1 - (1 - z) 2ct+1 ] (78) 

is such that Q(0) = and that 

^ = {2a + 1) [4*V(1 ~ z 2 y- x + ((1 - z) ff - (1 + zYf] > , Vz e [-1, +1], z Y 0, 

then we can conclude that dct M = if, and only if, z = and z is a zero of the 
function h(z). Then, as in the case of symmetric solution then ([77|) holds true. 
The Lemma is so proved. □ 

Remark 17. In fact, for symmetric solution a straightforward calculation gives 
that 

( 2- \ 




( E ' \ 




( P 


p' 







q' 




q 


a' 







V P J 




V 



(79) 



V J 

On the other hand, for asymmetrical solution corresponding to z = z as a straight- 
forward calculation gives also that 

/ Q(z) \ 
2"+ 1 h(z) > 
V2[(l+2) J -(l- Z )' T ] 2 (l-z 2 )VT^I 



/ E> \ 




/ 




1 \ 


P' 











q' 








-1 


a' 











V P 1 









/ 



2°+- i h(z) 

V2[(i+zr-(i-zr] 2 (i-z 2 )vr+z 

2°+ s h(z) 





V 



(80) 



where the function Q{z), defined in equation f7<S| ), is such that Q(—z) = —Q(z), 
Q(0) = and ^ > for any z £ (0, 1]. 

Now, we are ready to go back to the slope condition and to state the following. 

Lemma 9. There exists hs € (0, /12) such that for any h € (0, ^3) the following 
statements are satisfied. Let 

Fs(X) = UiJ 2 
where Y\ e is the symmetric stationary solutions. Then 



Tx F sW < 0. 
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Moreover, 

(i) Let a < (Jthreshold and let 

Fas(X) = UZW 2 

where ?/^ s e is the asymmetrical stationary solutions. Then 

lx Fas{x) < °" 

(ii) Let a > (Jthreshold and let 

F asl {X) = UTl\\ 2 and F as2 {X) = Ul s2 A 2 
where an d ^x* 2 are asymmetric stationary solutions. Then 

J-F osl (A)<0 and -^F as2 (A)>0. 
Proof. We consider, at first, the case of the symmetric stationary solution corre 

dE 
dij 



sponding to z x = z s = 0. In such a case from (1751) it follows that = 2 " > 

and thus dF ^^ < proving so the first statement. 

Now, we consider the case of asymmetrical stationary solution corresponding to 
z x 7^ 0. In such a case from (|80[) it follows that 

dE _ Q(z x ) 



dr) 2°+ 1 h{z x ) 

is an even function and where Q(z x ) ■ z x > 0. Hence, the sign of ^ only depends 
on the sign of h(z x ). We have then showed all the statements in Lemma[9l recalling 
that (see the results in Section |3|) 

If a < crthreshoid, then the asymmetrical stationary solution </>^ s e corre- 
sponding to z x > satisfies condition h(z x ) > 0; 
If a > a 'threshold, then the asymmetrical stationary solution 
sponding to z x > satisfies condition h(z x ) > 0; 

If (7 > crthreshoid, then the asymmetrical stationary solution <f>^ corre- 
sponding to z x > satisfies condition ft,(z A ) < 0. 

□ 

Remark 18. In the same way, the monotone decreasing behavior of 

FaW = U a xJ 2 

associated to the antisymmetric stationary solution follows. 

Finally, collecting the results of Proposition [T] and of Lemmata [5l [6] and [9] then 
Theorem [3] follows. 

Remark 19. In Theorem^ in case of a > crthreshoid and \r]\ — rj + , we did not 
obtain any conclusion about the orbital stability. Recall that rj + G (0,oo) is defined 
by r/ + = with z + € (0, 1) such that ri'(z + ) = (see Theorem^. Let 4>\+. e + 

be the corresponding asymmetric stationary solution to X + = f2 + uE + where 



E- 



v^W+* + [(— ) +(— ) 
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and e + is given by curj + jc. According to Remark\T^ we see formally 



dF(X) 



dX 



1 ■ .(i-W^ -1 



= 0, (81) 

\v\=v + 



a=a+ Cftd V drj J 

since rj / (z + ) = 0. Thus, we are required to prove the stability /instability for the case 
dF/dX — 0. In fact, this case would be included in (2) of Proposition^ and we 
would conclude that, when a > o 'threshold, at the transition point \n\ = rj + from 
(jj'xl to (f>1 S e ' we should have the instability. To show this fact exactly, it suffices 
to compute 4-p- and prove that it is not zero at X = X + , following the argument in 
Maeda [35] (see also some related conditions in O [H] )■ At least "formally" this 
may be seen as follows: we note that the use of Budan-Fourier Theorem ensures 
dX(z)/dz negative for h small. By formal calculations, 

d 2 F / W \V* r 1 /l \ .i_ 2 (d-q\ 2 1 ,i_ 1 d 2 r]} 

fx = <™% sH™! -*«>£}/(§)'• 

ua uz (ma v az az j \ az / 

We have seen in Section 3 that rf 1 (z + ) =/= 0, which implies 1 a=a+ 7^ 0. However 
a rigorous justification seems more complex and we do not pursue in this direction 
in the present paper. 

Appendix A. Stationary states for a non-linear toy model 

Here, we introduce, as a toy model, the semiclassical Schrodinger equation with 
two attractive symmetric Dirac's 5 which is partially investigated in |24) . 

ih^ = H ^ + eg\i>\ 2 ^, ||#,t)||=l, xeU,teR, (82) 

where 

H a = -h 2 -^+f38_ a + f3S +a 

for some a £ M. and f3 < 0. Hereafter, for the sake of dcfinitcness, we assume that 
,g=l. 

Even though this operator Ho with Dirac measures do not satisfy the assump- 
tions for the potential V(x) in the Introduction, the two- level approximation used 
in the previous sections is directly applicable to this example. In this section, we 
will give some remarks for the properties of Hq, and the general theory we have 
used in the previous sections, for example, Cauchy problem and the orbital sta- 
bility. We remark that a symmetric-breaking phenomenon for the cubic nonlinear 
Schrodinger equation with double Dirac potential is discussed in |21j too, but not 
in the semiclassical regime. 

A.l. Spectrum of the linear operator. The spectral problem 

d 2 



dx 

for (3 < is equivalent to the spectral problem 

H a ip = Eip 



ip = Sip 
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where we set E — E/ti 2 and where the linear operator 

d 2 

H a = ——^, with a = (3/h 2 , 
is self-adjoint on the domain 

D(H a ) = [il> G H 2 (R \ {±a}) n H\R) : ^-(±a + 0) - ?(±a - 0) = o^(±a) 
^ ax ax 

Let us recall some basic properties of the spectrum of H a (see, e.g., [H [53] for 
details.) 

The essential spectrum of H a is purely absolutely continuous and coincides with 
the positive real axis: 

o- e ss{H a ) = <7 &c {H a ) = [0, +oo) . 

The discrete spectrum consists of two eigenvalues, at least, given by means of 
the Lambert's special function W(x) such that W{x)e wtyX%> = x. 
If a < the discrete spectrum is not empty, in particular, 

- if a < — — , then the discrete spectrum of H a consists of only one eigenvalue 
Ei (a, a) defined as 

EAa, a) = — [W (-aae aa ) - aa} 2 ; 

- if a > — i, then the discrete spectrum of H a consists of two eigenvalues 
Ei(a,a) and £?2(a,a) where 

E 2 (a, a) = \r [W (+aae aa ) - aa} 2 . 

4a z 

The two associated eigenvectors take the form: 
i) Let 

k x = y/% = ^-[W {-aae aa ) - aa} 



e lklX x < — a 




2k\ -\-ia 
2k 

e +tklX X > +a 



t!S ( e - 4fc i x + e lklX ) -a<x<+a 



where C\ is the normalization constant given by 

Ifcil 



Ci 



ii) Let 



y/(2\k 1 \+a)(2\k 1 \a + aa + l) 



fE^=-L\W (+aae aa ) - aa} 
2a 
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where C 2 is the normalization constant given by 

Itol 



Co 



^-{2\k 2 \ + a) (2\k 2 \a + aa + l)' 
Remark 20. Recalling that the Lambert's special function W(x) has the following 



asymptotic behavior 



W{x) ~ x- x 2 + ^-x 3 + 0(x 4 ) 



2* 

then it follows that the splitting is exponentially small, namely 

\£ x - £ 2 \ ~ h 2 a 2 e aa = |! e <^A 2 = ^e-M" 1 '. 

Remark 21. The resolvent formula for H a is known: let h 6 Cq°(]R), k 2 € p(H a ), 
and $sk > 0. XTie resolvent is expressed as follows. 

({H a -k 2 ]- 1 h)(x)= f K a {x,y:k)h{y)d Vl 

JR 

with the kernel K a having the following form 

4 

K a [x,y;k) = K (x,y;k) + ^K 3 a {x,y]k) 



where 



K (x,y;k) = ^ klX ~ Vl 



«(2fc + id) ik\x+a\+ik\y+a\ 

2k((2k + ia) 2 + a 2 e l4ka ) 

_; rt 2 Jib 
ta K ik\x+a\+ik\y+a\ 



Ki(x,y;k) 

K °( X ' V ''Q 2k((2k + ia) 2 + a 2 e Ma )' 

^(z,^) = ^(-i.-y;*) 

K*(x,y;k) = Kli-x.-y-k) 



We consider here i/ie case a > — 1/a wit/i a < 0. In such a case we have that 
the linear problem has two negative non degenerate eigenvalues: 

Ex < E 2 < . (83) 

A. 2. Nonlinear problem. The local existence of solution in H 1 (M.), and conser- 
vation laws of energy and L 2 norm are verified in a similar way to [14] ; the authors 
in [T3] applied Theorem 3.7.1 of [5 to the case of a — 0. In our case, we take 
— H a + Ei for the operator A of Theorem 3.7.1 of [S]. Then this operator A is 
a self adjoint operator on X = L 2 (M.) with the domain D(A) = D(H a ), and also 
A < 0. We take Xa = H 1 (M.) whose norm is equivalent to norm 

IMIL = \\(d/dx)v\\ 2 + (1 - i^lMI 2 + a(\v(a)\ 2 + \v(-a)\ 2 ). 

Condition (3.7.2) of Theorem 3.7.1 of [5] is satisfied with p = 2, and other conditions 
hold since we are in one dimensional case. 

For the existence of bifurcation of stationary solutions, it suffices to repeat the 
similar arguments in Section 3 (Theorem [T]), but in i/ 1 (M) instead of H 2 (M.). 
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We can check the assumptions for the orbital stability/instability of stationary 
states in H (M), as in Section [5l using the two level approximation. However, due 
to the singularity of Dirac potentials, we cannot consider the linearized problem 
with a more smooth domain than H 1 (R) , as, for ex. , was considered in [12) . Remark 
also that H 2 regularity allows us simply to have the nonlinear instability assuming 
the existence of an unstable eigenvalue (e.g. [7]). We thus give some explanations 
here. 

We consider as follows the linearized problem around the real valued rescaled 
stationary state </> €; a (e and A are fixed here to discuss the general theory, so we 
denote it simply by <f> from now on) . 

dv 

— =Av + F(v), v = (ui,« 2 ) £ D(A) with v 1 =$lv, v 2 = 3w , (84) 
dt 

where A(vi, 1)2) — (L x ^ e V2, —L_i! e Vi). A is a linear operator in L 2 (R) with domain 

D(A) = fv £ H 2 (R \ {±a}) n H 1 (M) : -p(±a + 0) - -^(±a - 0) = aVj(±a), j = 

where H m (R) = H m (R) x H m (R) for m £ Z. The nonlinear term is given by 

F(v) = z{|0 + v\ 2a {cj) + v)- |0| 2ff+1 - (a + lM 2a v - a\<j>\ 2 °v}. 

This operator A generates its Co-semigroup on L 2 denoted by e tA . Concerning 
the spectrum of A, we have the following Lemma. We note that we complexify the 
space when we consider the spectrum problem of A. 

Lemma 10. a css (A) C iR. 

Proof. The operator A can be rewritten in the following form (still denoted by A 
with abuse of notation) 

Av = ~i{H a -\-(a + 1)|0| 2CT - <t|0| 2ct T} v 

where l~v — v is a non-symmetric bounded linear operator. We consider the 
operator iA as the operator Aq perturbed by the operator C, i.e. 

iA = A + C, 

where A = H a - A, and C = -{a + l)\<t)\ 2a - a\(j}\ 2a T. It suffices to prove that 
a ess (iA) C R. To this end, we remark the following facts. 

- Since 4> £ H l (M) C L°°(K), C is a bounded operator. 

- It is known that o~ ess (Ao) = [—A, +00) C R. 

- C [Aq + A + 1]~ is a compact operator; indeed, [Aq + A + 1]~ is an inte- 
gral operator with kernel given by K®(x, y; i) + Ylj=i Ka( x > V'i 0- *-*ne can 
see, for e.g., that K 3 a (x, y; i) (j = 1, 2, 3, 4) and \<j>\ 2a K°(x, y; i) are bounded 
on i 2 (R 2 , dx dy). This implies that C [A + A + l] -1 is Hilbert- Schmidt. 

Then a ess (Ao) — a ess (iA) by means of the Weyl criterion. □ 

As for eigenvalues of A, there are finitely many eigenvalues at the exterior of 
the essential spectrum for h small. Indeed, A < for h small. Our aim is now to 
conclude the following Proposition. 

Proposition 2. Assume that A has an eigenvalue A m with 3?A TO > 0, and that for 
any e > 0, there exists M > such that 

\\e tA v\k* < Me^&^'WvWw, (85) 
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for any v G L 2 (R) and for any t > 0. Then, there exists eq > 0, such that for any 
5 > there exist a time T and an initial data uq G D(H a ) satisfying \\uq — <^||jji < 5, 
and infflgR ||m(T) — e l6 cj)\\ > Eq. 

Proposition [2] means that the linearized instability implies the nonlinear insta- 
bility. 

We may prove Proposition [2] as in the proof of Theorem 6.1 of Part II of [19] 
or in 181. Note that we have the Dirac measures in the equation and we do not 
expect that the solution is smooth as we have mentioned before, thus we make use 
rather of the time derivative, mimicking the proof of [8], than of the way of [19]. 
Here, for the sake of completeness, we give an outline of proof. 

Proof. (Sketch of proof) Let z m be the associated eigenfunction to A m . Let us(t) 
be the solution of ([521 with initial data u$(0) = 4> + 5z m . Since 0, z m e D(H a ), 
u s (-) £ C([0,T],D(H a )) n C^lO^lL 2 ) for some T > (see Theorem 3.1 of [T]). 
Remark that ug(t) — e~ lXt ((f> + vs(t)) with v$(t) satisfying ([84]) with vs(0) — Sz m . 
vs(t) satisfies the following integral equations for any t g [0,T], 

v B (t) = 6e x '"z m + f ' e iT ^ A F(v s (T))dT, 
Jo 

d t v s (t) = \ m 6e x ™z m + e tA F(6z m )+ f e^ A d T F(v 5 (T))dT. 

Jo 

Since we are in the one dimensional case, it is easy to estimate the nonlinear term 
F(v) as follows. 

11^(^(^)11 < C(\\v s (t)\\ Hl +\\v s (t)\\%)\\d t v s (t)l 
\\F(v s (t))\\ < C 2 (\\v s (t)\\ H i + Wvsm 2 ^ 1 ). 

Then, for some Co > and for some Tg when 5 is sufficiently small, we may estimate 

\\v S (t)\\„i+\\d t v s (t)\\ < 2C <5e A '"*, 

for any t € [0, Ts] . We apply this quantity || vs (t) ||#i + 1| d t vs (t)\\ as Vs (t) in Theorem 
2 of [8j. We then repeat their arguments in [8J to get ||u<s(T5)|| > (<5/2)||z m ||. □ 

We complete our whole arguments with a verification of the existence of an 
eigenvalue A m satisfying ([85]) . It follows from jT8] or Part I of [19] that there exists 
a non-zero real eigenvalue of the linearized operator A, if (2) or (3) of Proposition 
[T] in Section [5] hold. Let Ao be the maximal positive eigenvalue. Once we have 
proved the spectral mapping theorem a(e At ) — e"^ 1 , the spectral radius of e At is 
gA t Thus we have ([85]) using Lemma 3 of [30] ■ This implies that we can take Ao 
as X m in Proposition [2] 

The spectral mapping theorem in fact follows from a resolvent estimate in Lemma 
llll below. combined with the arguments in |16j . 

Lemma 11. Let z = a + it with a, r € K and a ^ 0. For |r| sufficiently large, 
there exists a constant C a > 0, such that 

\\(z - A)- 1 ^) < C a . 
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Proof of Lemma XTlS. We begin with some preparations. For fixed z = a + it with 
a £ R \ {0} and r £ R, we write the operator z — A as follows. 



z -H a 
H a z 



<p 2a + A 

-A - (2cr + f )</> 2<T 



= M^Id-M^BxA 



Indeed, we see that z iR, therefore, by Remark [211 the inverse of H 2 + z 2 = 
— iz)(H a + iz) exists, thus the inverse of M z exists too. We can express M" 1 
as follows. 

M -i_( z{(H a ) 2 + z 2 }^ H a {{H a f + z 2 }-^ 
Mz ~ \-H a {{H a f + z 2 }' 1 z{(H a ) 2 + z 2 }-\ 

We estimate now the inverse M~ x by means of the following Lemma. 

Lemma 12. Let a ^ 0. There exist C a , tq > such that for any z = a + it with 
\t\ > tq, we have 

C a 



m: 



|£(L2) 



< 



1+ T 



Proof of Lemma \12{ We benefit from the explicit resolvent formula of H a in Re- 
mark [2D Let 

f a (x) = {[H a -k 2 ]- 1 h){x) 



and consider k 2 = iz = — r 



First, we remark that 



/„(*) = J K (x, y; k)h{y)dy = ± (e* fc H * h(-)) (x) 
may be estimated, by Young inequality, as follows. 



1/0,1 ~2\k\ 



,ik\-\ 



*M0 



< 



2|fe| 



ik\-\ 



TA / 1 



+ + Ot~ 3/2 , as |r| -> oo , 



and 3Vt + za ~ |^ for r > 1. 



P a (x)= / Ki(x,y;k)h(y)dy, j = 1,2,3,4. 



Next, we set 



By this definition, / Q = /o + J2j=i fa- Thus we estimate each term f a . For 
example, 

a(2k + ia) 



/«(*) = 
, for suffi 

im < 



2k((2k + ia) 2 + a 2 e lika ) 
and then, for sufficiently large |r|, 
C 



ifc|a;+o| / ik\y+a\ 



h(y)dy 



< 



\k\ 2 \Sk\ 
C 



Jk\y+a\ 



ik\-+a\ 



h(y)dy 

N < L 



< 



c 



|/fc| 2 |3fc| 



L 1 



c 
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Similarly, the other terms /£, j — 2,3,4, are estimated. Thus, it follows that for 
3z = a fixed and SRz = — t large enough, then 

\t\ 

since H a is a self-adjoint operator. Therefore, decomposing H a {(H a ) 2 + z 2 }~ 1 as 

H a {{H a f + z 2 }- 1 = {H + iz)- 1 + iz{H - iz)-\H + iz)' 1 , 
we also obtain, for large |t| ^> 1, 



l^ a {(jj a ) 2 + ^ 2 }- 1 IU(^)< T ^- 

1+1 



Similarly, for large |r| 



||4(iJ a )2 + ^ } -i| U(L2) <_^ 



□ 



We go back to the proof of Lemma [TT] We put T z = M z 1 B E: \, and we write 
entries of this operator T z : 

_ (H a {(H a ) 2 + z 2 }-\-\ - (2a + 1)0 2 -), z{(H a ) 2 + z 2 }~\cf> 2 * + A) 
lz \z{{H a f + z 2 }- l {-\-{2a + l)^ a ), -H a {(H a ) 2 + z 2 }-\^ + \) 

bmce we are m one dimension, it follows that 4> e H 1 ^) C L°°(R), thus we can 
estimate, for example, as 

\\H a {(H a ) 2 + z 2 }- 1 ^ + X)\\ £m < C\\H a {{H a f + z 2 }- 1 ^^. 

Therefore, combining with the above proof for Lemma 111) we have that for any r 
with |t| > To, ||2«||£(l 2 ) < 1/2- This implies immediately for any u G L 2 

\\{Id-T z )u\y > \\u\y - \\T x u\\u > (l/2)||u|| L 2, 

that is, Id—T z is invertible for |r| > To. Then, finally, we get that for any z = a+ir 
with |t| > To, a ^ 0, 

Uz-A)- 1 ^ = Uld-T^M-'Wc^ 

< \\(Id-T z )- 1 \\ c{h2) \\M z 1 \\ c{h2) <2C a . 

The proof of Lemma fTTI is then completed. □ 

Lastly, recall that the assumptions (2) or (3) of Proposition[T]in Section [5] ensure 
the existence of a positive real eigenvalue of A. As we checked in Section [5J the 
assumptions (2) or (3) of Proposition [1] in Section [5] may be verified, for small 
h > 0, depending on a, 77, and the sort of stationary solution. Namely, Theorem [3] 
in Section [5J is valid for Eq.(|82]). 
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